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ABSTRACT 


A  generalized  automated  data  acquisition  system  was  designed  for  the  Naval  Postgraduate 
Set  ol  Low  Speed  Wind  Tunnel.  A  specific  application  of  this  system  was  to  upgrade  the  current 
"Cylinder  Drag  Experiment”  conducted  during  AA2801  "Aero  Laboratories  I“,  an  introductory 
aeronautical  laboratory  course  taught  at  the  Naval  Postgraduate  School.  Two  methods  of  drag 
determination  were  used:  pressure  distribution  and  wake  analysis  (momentum  method).  Data  from 
these  two  methods  were  collected  by  a  system  based  on  a  high  speed  analog/digital  computer  board, 
a  standard  486  IBM-type  PC  and  data  acquisition  software.  Characteristic  methods  of  reducing  data 
from  this  experiment  are  discussed.  The  results  obtained  by  analyzing  the  acquired  data  compared 
favorably  to  empirical  data  from  previous  circular  cylinder  coefficient  of  drag  experiments.  This 
automated  data  acquisition  system  will  facilitate  future  research  and  instructional  use  of  the  wind 
tunnel . 
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TABLE  OF  SYMBOLS  AND  ABBREVIATIONS 


p0  =  total  (stagnation)  pressure 

p  =  static  pressure 

q  =  dynamic  pressure  =  p0  -  p 
pm  =  freestream  density 

(po)*  =  freestream  total  pressure  (associated  with  freestream  velocity) 
p*  =  freestream  static  pressure  (associated  with  freestream  velocity) 
q„,  =  freestream  dynamic  pressure  =  (po)„  -  p. 

(Po)w  =  wake  total  pressure  measured  behind  the  cylinder  at  the  probe 
((Po)w  assumed  to  =  (po).  when  taken  outside  wake) 
pw  =  static  pressure  measured  behind  the  cylinder  at  the  probe 
(pw  assumed  to  =  p„.  when  taken  outside  wake) 
qw  =  wake  dynamic  pressure  =  (p0)w  -  pw 
pcyl  =  pressure  measured  at  the  port  on  the  cylinder 
Cp  =  pressure  coefficient 

CD  =  drag  coefficient  (3-D) 

cd  =  drag  coefficient  relating  to  total  drag  or  profile  drag 

cdp  =  drag  coefficient  elating  to  form  drag  only 

cdu  =  uncorrected  drag  coefficient 
6  =  azimuth  of  the  rotator  (signal  from  potentiometer) 

y  =  position  of  the  traverse  (signal  from  potentiometer) 
r  =  radius  of  cylinder 

R  =  gas  constant 

d  =  diameter  of  the  cylinder  =  2  x  R 

da  =  incremental  area 

df  =  incremental  pressure  drag  force 

dd  —  incremental  form  drag  component 

dl  -  incremental  induced  drag  component 

d  =  drag  (2-D) 

1  =  length 

w  =  width 

h  =  height 

C  =  cross-sectional  area 
F  =  measured  force 

c  =  reference  chord 

S  =  reference  area 

A  =  area  (2-D) 

V  =  volume  (3-D) 

e  —  blockage  factor 

V„  =  freestream  velocity 

Vw  =  velocity  as  measured  along  the  traverse  behind  the  cylinder 
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vmu>  =  the  velocity  into  a  closed  system 

V^,  =  the  velocity  out  of  a  closed  system 

D  =  total  drag  unless  specifically  defined  as  "profile  drag" 

D(  =  shear  drag  -  drag  due  to  skin  friction 

Dp  =  form  drag  -  drag  due  to  separation 

Re  =  Reynolds  number 

eff  =  (subscript)  effective  (corrected  value) 

u  =  (subscript)  uncorrected  value 


Conversion  Constants 

latitude  correction  =  -0.0245  inches  Hg 

temperature  correction  (see  Table  4-2  on  bulletin  board  near  tunnel) 
=  3.8  X  10 7  (lbf-sVft2 
R,ir  =  53.3  (ft-lbf)/(lbm-°R) 

°R  =  °F  +  460 
1  ft  =  30.48  cm 
1  cm  H20  =  0.0142  psi 
1  psi  =  2.036  inches  Hg 

Qicst  icction  =  (Ap  *  0.243)/. 895  (NPS  wind  tunnel  calibration  equation) 
TF  =  1.04  (NPS  wind  tunnel  turbulence  factor) 
gc  =  32.174  ft/ sec2  (universal  gravitational  constant) 
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I.  INTRODUCTION 


A  classic  experiment,  drag  measurement  over  a  circular  cylinder  in  a  wind  tunnel, 
serves  as  an  excellent  primer  on  wind  tunnel  use,  data  acquisition  and  basic 
aerodynamics.  This  thesis  served  to  upgrade  an  existing  AA2801  (Aero  Laboratories  I) 
experimental  lab  (see  Figure  1)  while  further  exploring  the  theory  behind  the  methods 
used  to  measure  drag.  Chapter  n  examines  the  Naval  Postgraduate  School’s  Low  Speed 
Wind  Tunnel,  including  a  detailed  procedure  for  calibrating  the  wind  tunnel.  Chapter 
m  will  discuss  the  theory  behind  drag  measurements  around  a  circular  cylinder.  Chapter 
IV  compares  the  former  and  new  procedures  used  in  the  AA2801  laboratory.  The  results 
of  the  new  experiment  using  actual  class  data  are  then  analyzed  in  Chapter  V  and  some 
conclusions  are  drawn  in  Chapter  VI.  Appendix  A  contains  the  former  AA2801  Cylinder 
Drag  Experiment  Handout.  Appendix  B  contains  a  copy  of  the  upgraded  handout  that 
was  used  on  9  November  1993  with  a  class  of  graduate  students.  Appendix  C  contains 
a  readout  of  the  experimental  data  obtained.  Appendix  D  is  the  MATLAB  script  file 
"DRAGl.m"  used  to  analyze  data  from  the  upgraded  lab.  Appendix  E  has  several 
MATLAB  script  files  for  creating  velocity  tables  and  curves  which  were  derived  from 
wind  tunnel  calibration  data. 


1 


Figure  1.  The  Test  Section  Showing  The  Wake  Traverse  and  Probes  Apparatus 
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n.  N.P.S.  LOW  SPEED  WIND  TUNNEL 


A.  DESCRIPTION  AND  OPERATING  PROCEDURES 

The  NPS  Laboratory  Manual  for  Low  Speed  Wind  Tunnel  Testing  [Ref.  1] 
provides  an  adequate  discussion  of  the  basic  system  description  of  the  NPS  Low  Speed 
Wind  Tunnel  and  its  chief  operating  procedures  and  should  be  read  thoroughly  by  all 
users.  What  follows  is  a  summary  of  this  description  and  those  procedures  of  interest 
to  the  Cylinder  Drag  Experiment. 

Figure  2  depicts  the  layout  of  the  NPS  Low  Speed  Wind  Tunnel  and  the  location 
of  the  test  section.  The  wind  tunnel  is  a  closed-circuit,  single-return  type.  It  is  powered 
by  a  100  hp  electric  motor  which  drives  a  three-blade  variable  pitch  fan. 

A  four-speed  common  truck  transmission  is  used  in-line  with  the  motor  and  allows 
the  wind  tunnel  to  achieve  tunnel  speeds  up  to  200  MPH.  The  Cylinder  Drag 
Experiment  in  AA2801  operates  at  approximately  1 15  MPH  and  for  this,  fourth  gear  is 
used.  The  contraction  cone  prior  to  the  test  section  has  a  10:  l  contraction  ratio  which 
serves  to  increase  the  velocity  of  the  airflow  as  well  as  causing  the  flow  to  become  more 
uniform  entering  the  test  section. 
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Figure  2.  NPS  Low  Speed  Wind  Tunnel 


The  dimensions  of  the  test  section  are  given  in  Table  1 . 
TABLE  1.  Test  Section  Dimensions. 


length  (1)  . 48.0  inches 

width  (w) . 45.0  inches 

height  (h) . 32.0  indies 

total  cross  sectional  area  (C)  . 9.9  sq.  ft 

above  the  reflection  plate: 

height  (h^n) . 28.4  inches 

cross  sectional  area  ( C . 8.8  sq.  ft 


Figures  3  and  4  are  views  from  inside  the  test  section.  The  wind  tunnel  uses  a 
water  micromanometer  to  reference  tunnel  velocity  by  measuring  Ap,  the  difference 
between  two  static  pressures:  p,  at  the  beginning  of  the  contraction  cone  and  p2  at  the  exit 
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Figure  3.  Test  Section  View  -  Into  the  Freestream 


of  the  contraction  cone.  Since  the  airspeed  is  increased  between  the  beginning  and  exit 
of  the  cone,  there  is  a  decrease  in  the  static  pressure  as  the  air  seeks  to  return  to  the 
lower  ambient  static  pressure.  This  equalization  is  due  to  a  breather  slot,  which  allows 
ambient  air  to  enter  the  circuit  and  mix  with  air  leaving  the  test  section.  Hence,  Ap 
changes  primarily  with  increasing  (or  decreasing)  p, ,  as  p2  remains  essentially  constant. 
Figure  5  shows  the  water  micromanometer.  Through  calibration  and  associated  equations 
and  charts  (see  Section  C),  this  Ap  is  related  to  airflow  velocity. 
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Table  2  illustrates  the  operating  instructions  used  during  this  lab. 


TABLE  2.  Wind  Tunnel  Operating  Procedures 


STEP 

PROCEDURE 

1 

Decide  on  a  required  tunnel  speed;  convert  to  Ap  (H20  cm) 

2 

After  confirming  micromanometer  reads  0.00  with  tunnel 
fully  stopped,  set  selected  Ap  in  window 

3 

Check  tunnel  doors  closed  and  test  section  free  from  FOD; 
check  transmission  for  proper  gear;  turn  power  switch  to  ON 

position 

4 

Depress  and  hold  START  button  until  green  light  goes  out 
momentarily  (meaning  fan  is  turning  and  set  to  minimum 
pitch);  release  START  button 

5 

Wait  5-10  seconds;  depress  COARSE  UP  button  until  water 
rises  to  enter  the  crosshair  ring  on  micromanometer;  toggle 
FINE  UP/DOWN  to  maneuver  meniscus  to  rest  on  crosshair; 

tunnel  is  set 

STOP 

Toggle  COARSE  DOWN  until  water  level  is  below  the  10  cm 
mark  then  momentarily  depress  STOP  (button  with  red  light); 
secure  power  switch;  enter  time  into  tunnel  logbook 

EMERGENCY  STOP: 

Secure  switch  on  red  box  near  breather  slot 

B.  CORRECTION  FACTORS 

Data  taken  from  within  the  wind  tunnel  must  first  be  corrected  and  calibrated  due 
to  the  physical  geometry  of  this  particular  wind  tunnel.  Of  primary  importance  are  the 
actual  Reynolds  number  and  dynamic  pressure  being  experienced  by  an  object  inside  the 
test  section.  Since  these  are  functions  of  velocity,  an  accurate  method  of  obtaining  test 
section  velocity  must  be  available.  This  velocity  may  be  derived  from  the  difference  p,  - 
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Figure  5.  N.P.S.  Low  Speed  Wind  Tunnel  Water  Manometer 
p2,  but  as  was  mentioned  earlier,  the  tunnel  Ap  is  meaningless  unless  it  has  first  been 
calibrated  to  a  known  standard.  This  calibration  procedure  (which  should  be  repeated 
before  any  important  research  is  undertaken),  is  simply  the  tunnel  calibration.  The  actual 
procedure  and  results  are  discussed  in  the  next  section.  Here  it  is  pointed  out  that 
Ap,once  properly  calibrated,  may  then  be  used  to  provide  the  test  section  freestream* 

’"Freestream"  will  henceforth  refer  to  any  parameter  measured  in  the  test  section 
which  is  unaffected  by  any  object  (eg.  an  airfoil)  placed  in  or  any  mechanism  directed 
into  (eg.  a  jet  of  air)  the  test  section,  and  will  be  shown  with  (<»)  as  the  subscript. 
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velocity  V,,  as  demonstrated  by  equation  (1).  The  development  of  the  tunnel 
calibration.equation  is  explained  in  the  next  section. 


A  p  =  ( calibration  factor)  x  q  =  (calibration  factor)—  pmVj 

2 


(1) 


Besides  calibrating  for  q,  there  is  a  correction  for  the  turbulence  within  the  test 
section  called  the  "Turbulence  Factor,"  or  TF.  It  is  found  experimentally  and  was  not 
repeated  for  this  thesis  [Ref.  l:pp.  32-33]. 


385,000 

^ eCUT 


(2) 


The  TF  was  found  experimentally  for  the  NPS  wind  tunnel  to  be  1.04.  The  TF 
effectively  raises  the  Reynolds  number. 

The  final  corrections  of  interest  to  this  lab  are  the  tunnel  boundary  corrections  due 
to  solid  and  wake  blockage  caused  by  the  presence  of  the  cylinder  in  the  test  section. 
These  corrections  account  for  the  test  section  walls  and  the  apparent  "closed  system"  the 
data  acquisition  system  experiences.  These  are  corrections  for  two-dimensional  testing. 
Solid  blockage  correction  a*  compensates  for  increased  velocity  and  dynamic  pressure 
due  to  the  blockage  caused  by  the  model  [Ref.  l:p.  36]. 
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= 


K ^Volume 


3 

"»  2 

'<ar 


(3) 


K,  is  0.S2  for  a  model  spanning  the  tunnel  height  and  Cefr  is  the  effective  cross-sectional 
area  of  the  test  section. 

The  second  factor  is  £wb,  or  the  wake  blockage  correction.  The  decrease  in  velocity 
caused  by  a  wake  results  in  an  increased  velocity  outside  the  wake  as  the  "closed  system" 
tries  to  maintain  a  constant  mass  flow  rate  within  the  test  section.  This  results  in  an 
increased  dynamic  pressure  on  the  model.  [Ref.  l:pp.  36-37] 


*wb 


(4) 


Here,  d  is  the  diameter  of  the  cylinder,  w  is  the  width  of  the  test  section  and  C*  is  the 
uncorrected  drag  coefficient.  Reference  1  (pp.  36-37)  goes  into  broader  detail  about 
what  Cdu  represents.  For  the  purposes  of  this  experiment,  C*.  will  be  chosen  to  be  c„ 
(and  thus  labeled  "c^")  found  by  each  of  two  methods  described  in  the  next  chapter. 
Reference  1  offers  a  simple  estimate  for  total  blockage  on  page  37,  so  no  significant  loss 
in  precision  will  be  experienced  in  calculating  cd. 


Cd  ~  Cd(^  ^esb  ~  ^Ch 


(5) 
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Thus  [Ref.  l:p.  37], 


K  ^Volume  4 

ttoUi  =  3  +  ~^Cd  (6) 

C  2 

£t0Ui  is  then  used  in  further  corrections  to  velocity  and  the  Reynolds  number  [Ref.  1  :p. 
38], 


V  =  K„(l  + 


(7) 


R^eff  Re  +  E  »*•()"(  7!F) 


(8) 


C.  CALIBRATION  AND  RESULTS 

The  relationship  between  Ap  and  tunnel  q  is  essentially  linear  and  provides  the 
tunnel  calibration  curve.  The  slope  and  intercept  of  the  calibration  curve  comprise  the 
tunnel  calibration  factor  which  is  the  object  of  the  calibration. 

The  tunnel  was  calibrated  on  26  April  1993  using  three  steps: 

1.  The  output  from  a  digital  manometer  was  compared  to  that  of  a  known 
standard  and  a  correction  factor  was  determined  for  the  digital  manometer. 

2.  The  reading  from  the  water  micromanometer  used  for  the  wind  tunnel  (Ap) 
was  compared  to  that  of  the  digital  manometer,  which  now  read  wind  tunnel 
q  using  a  portable  pitot-static  tube  placed  in  the  center  of  the  test  section. 

3.  The  data  obtained  were  reduced  to  a  plot  of  Ap  versus  corrected  tunnel  q; 
this  plot  was  the  calibration  curve  sought. 
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1.  Calibration  Step  One 


The  digital  manometer  was  simply  a  transducer  with  one  pressure  port  labeled 
"high”  and  one  pressure  port  labeled  "low."  The  high  port  was  connected  to  the  output 
of  a  known  standard,  in  this  case  a  device  called  the  "Schmidter. M  Next,  the  manometer 
reading  of  the  Schmidter  was  compared  to  the  output  of  the  digital  manometer  over  the 
pressure  range  necessary  for  the  test.  The  "Schmidter"  is  shown  in  Figure  6,  and 
includes  a  cylinder  which  can  be  manually  pressurized  and  a  water  manometer  reading 
in  centimeters. 


Two  tests  were  made  and  the  results  were  then  averaged  and  plotted.  Figure  7  shows 
the  resulting  conversion  derived  from  this  comparison.  The  subsequent  slope  and 
intercept  (offset)  of  the  plotted  line  were  then  applied  to  data  found  during  step  two. 
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Three  things  are  to  be  noted  concerning  Figure  7:  (1)  the  digital  manometer  was 

designed  to  read  out  in  inches,  not  centimeters,  and  thus  the  conversion  shown  in 
Figure  7  shows  inches;  (2)  it  can  be  seen  from  the  calibration  equation  in  Figure  7  that 
the  digital  manometer  had  a  linear  error.  This  explains  why  the  digital  manometer  had 
to  First  be  calibrated  to  a  known  standard  before  it  could  subsequently  be  used  to 
calibrate  Ap;  and  (3)  finally,  since  the  digital  manometer  used  had  an  operating  limit  near 
13  centimeters,  higher  runs  were  not  pursued.  Even  though  the  calibration  was  not 
ronducted  at  higher  pressures,  other  data  are  available  w  hich  ensure  that  the  calibration 
factor  remains  linear  for  Ap’s  above  13  centimeters.  The  AA2801  Cylinder  Drag 
Experiment  uses  the  wind  tunnel  in  4th  gear  at  15  cm  HzO  Ap. 

2.  Calibration  Step  Two 

The  next  step  involved  placing  a  standard  pitot-static  tube  into  the  wind 
tunnel,  extending  approximately  10.5  inches  from  the  top  center  of  the  test  section  and 
facing  into  the  freestream.  The  tunnel  was  operated  in  increments  of  1  cm  HjO  (tunnel 
speed  gauge  (Ap))  and  air  velocity  (in  MPH)  was  recorded  from  the  tunnel  airspeed 
indicator.  At  each  increment,  the  digital  manometer  output  of  the  pitot-static 
tube  was  recorded  (see  Table  3).  The  tunnel  was  operated  first  in  4th  gear,  and  then  the 
operation  was  repeated  in  3rd  gear.  The  values  read  by  the  digital  manometer  were  then 
averaged.  The  correction  from  Figure  7  was  applied  to  this  average  and  then  plotted  in 
Figure  8. 
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Figure  7.  Data  From  Step  One  of  the  Calibration 
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TABLE  3.  Step  Two  Results 


First  run  in  4th  gear/Second  run  in 
3rd  gear 

i 

Corrected  q  (cm) 

AP 

(cm  H20) 

Average  Digital 
Manometer 
(inches) 

1.00 

0.41 

39 

0.88 

2.00 

0.87 

47 

2.01 

3.00 

1.33 

56 

3.13 

4.00 

1.74 

63 

4.13 

5.00 

2.20 

69 

5.25 

6.00 

2.68 

75 

6.42 

7.00 

3.13 

79 

7.52 

8.00 

3.59 

85 

8.64 

9.00 

4.06 

90 

9.78 

10.00 

4.51 

95 

10.88 

11.00 

5.02 

99 

12.12 

12.00 

5.44 

103 

13.15 

3.  Calibration  Step  Three 

Step  three  was  to  plot  the  first  and  fourth  columns  of  data  from  Table  4 
against  each  other.  Again,  a  simple  curve  fitting  program  was  used  to  find  the 
calibration  curve.  Here,  the  slope  and  offset  of  Figure  8  provide  the  tunnel  calibration 
factor,  and  the  equation  is  0.243  +  0.895  *  q,  with  q  given  in  cm.  of  water.  This 
differs  from  an  earlier  calibration  of  0.93  *  q.  It  should  be  noted  that  the  two  formulas 
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are  merely  curve  fits  to  measured  data,  and  do  not  vary  as  much  as  might  be  supposed. 
The  curves  produce  identical  values  at  about  7  cm.,  and  at  IS  cm.  only  differ  by  about 
0.25  cm.  At  10  cm.  the  difference  is  about  0. 1  cm.,  well  within  the  fluctuation  level  of 
the  micromanometer  reading  during  the  periodic  pressure  surges  in  the  tunnel. 

NPS  Tunnel  Calibration  Factor  for  the  NPS  Aerolab 
Low  Spead  Wind  Tunnel  Howard/Millar  - 
Tunnel  'Delta  4/2«/t3...plot  revlaad  12/3/93) 

P*  (cm  H20) 

13 
12 
11 
10 
9 
6 
7 
6 
5 
4 
3 
2 
1 
0 


Figure  8.  Calibration  Curve  from  Step  Two  data 
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Dynamic  pressure  for  incompressible  flow  is 


tunnel  q  =  qm 


(9) 


and  Ap  (in  lb/ft2)  is 


Ap  =  [p,  (cm)  -  p2  (cm) 


(10) 


where 


*  62.35 


ft3 


(11) 


The  wind  tunnel  calibration  equation,  from  Figure  8,  is 

Ap  -  0.243  (cm  Hfi)  ♦  0.895  (12) 


And  combining  equations  (9)  through  (12)  and  taking  density  in  slugs/ft3. 


i - 

II 

/ — \ 

u* 

2(Ap  -  0.243) 

0.895p„ 

This  equation  has  been  written  into  a  MATLAB  code  which  uses  different  unit 
conversion  factors  to  create  updated  velocity  tables.  This  code  is  contained  in  Appendix 
E. 
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ffl.  CYLINDER  DRAG  THEORY 


A.  AERODYNAMIC  DRAG  CLASSIFICATION 
From  Reference  2: 


TABLE  4.  Drag  Classification. 


TOTAL  DRAG  -  D  j 

SHEAR  -  Df 
(due  to  skin  friction) 

Dominates  streamlined 
bodies 

NORMAL 
(due  to  pressure) 

Form  -  Dp 
(due  to  separation) 

Dominates  blunt  bodies 

Induced 
(due  to  lift) 

Profile  Drag  =  Df  +  Dp  (2-D  flow) 

Parasite  Drag  =  Df  +  Dp  (3-D  flow) 

Though  not  used  in  this  experiment,  a  wind  tunnel  balance  may  be  used  to  directly 
measure  the  drag  (and  lift/side)  forces  on  a  cylinder  in  the  test  section.  The  balances 
used  are  divided  into  two  categories:  internal  (sting-type)  and  external.  The  NPS  Low 
Speed  Wind  Tunnel  has  provisions  for  an  external  balance  using  a  strain-gage  bridge 
network  which  measures  lift  and  drag. 

Methods  which  use  pressure  to  find  drag  include  the  pressure  distribution  method, 
which  finds  form  drag  only,  and  wake  analysis,  which  finds  profile  drag. 


B.  MEASURING  DRAG  BY  PRESSURE  DISTRIBUTION 
1.  Theory. 


Figure  9.  Pressure  Forces  on  the  Surface  of  a  Cylinder 

As  shown  in  Figure  9,  which  was  derived  from  a  similar  figure  in  Reference  3  (p. 
58),  there  is  a  pressure  force  df  acting  normal  to  the  surface  of  an  airfoil.  There  is 
another  pressure  force  not  shown  which  acts  tangentially  to  the  surface.  This  tangential 
force  is  called  shear  and  is  a  result  of  skin  friction.  However,  for  blunt  bodies,  drag 
related  to  surface  pressure  dominates  and  the  tangential  pressure  force  is  commonly 
ignored  [Ref.  3:p.  60].  So,  the  force  shown  is  the  significant  contributor  and  has  a 
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component  called  form  drag  dd  associated  with  the  streamwise  drag  acting  on  the 
cylinder.  This  form  drag  may  be  found  by  starting  with  the  normal  pressure  force  as 
follows: 

df  =  p  da  =  p  (r  dd)  (H) 

where  df  is  acting  over  a  unit  span  and  r  is  the  radius  of  the  cylinder,  as  shown.  It 
follows  that 

dd  =  df-cosd  =  prd&  cos0  (15) 

where  dd  is  the  form  drag  component  and  dl,  shown  in  Figure  9,  is  the  induced  drag 
component.  By  integrating  over  the  whole  surface  of  the  cylinder  to  sum  all  the  dd 
components, 

f  2*dd  =  rf2np  cos  Qdd  (16) 

Jo  Jo 


which  equals  the  form  drag ,  or 


Dp  =  rj2np  cosddd 


(17) 
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The  static  pressure  p  acting  normal  to  the  cylinder’s  surface  can  be 
considered  as  a  local  static  pressure  which  varies  with  the  pressure  port  location  (as  the 
cylinder  is  rotated  from  6  =  -20°  to  6  =  +180°)  and  the  freestream  static  pressure 
which  is  constant.  The  freestream  static  pressure  p„  can  be  pictured  as  a  "blanket"  of 
constant  thickness  enveloping  all  360°  of  the  cylinder.  Integration  of  this  "blanket"  over 
a  closed  surface  will  result  in  zero. 

As  noted  previously,  integration  of  the  local  static  pressure  over  the  cylinder  acting 
in  the  downstream  direction  will  produce  the  form  drag.  Calling  this  local  varying  static 
pressure  pcyl, 

Dp  =  r  f**Pcypos&dQ  (18) 

And  due  to  symmetry, 

[2Kcosddd  =  0  (19) 

Jo 

which  can  be  written  due  to  constant  p„, 

rj2*pjcos&dd  =  0  (20) 

Therefore,  the  part  of  the  pressure  which  contributes  nothing  to  the  form  drag  can 
be  removed  from  the  previous  expression: 

Dp  =  rj2n  Pcyl  cos 0d0  -  rj'2n  Pm  cos QdQ  (21) 


21 


which  simplifies  to 

Dp  =  rfj9  (p^  -  pm)  cos  Odd  (22) 


Furthermore,  since  the  drag  equation  is  already  per  unit  span,  it  can  now  be  made 
dimensionless.  Also,  as  two-dimensional  flow  is  being  examined,  the  total  reference  area 
S  reduces  to  the  chord  length  c  x  1 ,  which  for  the  cylinder  is  2r  x  1  and 


C-,= 


r/02*(Pcy  "  PJ  cosBdB 


qm  2 r 


(23) 


The  r’s  cancel,  leaving 


^f0UKt  -  P->  cosQdB 


(24) 


The  pressure  coefficient,  Cp,  is  a  dimensionless  ratio  between  the  local  pressure  acting 
normal  to  a  given  surface  minus  the  freestream  static  pressure  (p.),  and  the  freestream 
dynamic  pressure  (q*).  With  reference  to  equation  (24),  since  q„  at  a  fixed  tunnel  speed 
is  simply  a  constant, 

c  =  lf2n(  ZsL Hz  )  cosed©  (25) 

*9  2Jo  qm 
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which  is  the  same  as 


c,  =  lf2*C  cos QdB  (26) 

dr  2J0  * 


The  cosine  function  is  an  even  function  and  by  symmetry  the  last  equation  can  be  written 
as 


The  following  authors  arrive  at  the  same  conclusions  using  slightly  differing 

approaches: 

1.  Anderson,  [Ref.  2:p.  209,  equation  3.130] 

2.  Bertin  and  Smith,  [Ref.  3:p.  59,  equations  2.40,  2.42] 

3.  Rae  and  Pope,  [Ref.  4:  p.221,  equation  4.29] 

4.  Hoemer,  [Ref.  5:p.  1-9,  equations  4-6] 

There  may  be  some  concern  regarding  the  sign  convention  used  in  this 
experiment  as  presented  here.  In  the  AA2801  Cylinder  Drag  Lab,  it  is  assumed  that  8 
=  0  and  8  —  tt  are  located  as  shown  in  Figure  9.  Thus,  integrating  an  otherwise 
dominating  positive  area  yields  a  correspondingly  positive  c^.  Classic  discussion 
normally  has  the  limits  of  integration  reversed,  with  8  =  0  beginning  at  the  leeward  side 
of  the  cylinder  and  growing  positive  in  a  counterclockwise  fashion.  Additionally,  after 
seeing  a  plot  such  as  Figure  10,  which  compares  theoretical  (inviscid)  Cp  to  actual 
(viscous)  Cp,  one  might  be  led  to  believe  that  a  negative  sign  would  be  required  in  front 
of  Equation  (27)  just  to  make  c#  positive. 
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Figure  10.  Inviscid  C'r  ns  Compared  to  Actual  C,, 

Actually,  the  integrand  is  not  simply  C,.:  it  is  Cr  v  cos 0  Figme  I  I  makes  this  appaient 
and  addresses  the  positive  area/positive  sign  question 


Figure  II.  Cr*cosfl  compared  to  C,. 

(  .  MKASl'RIN(»  I)RA(;  BY  WAKK  ANALYSIS  (MOMLN 1  I'M) 

1 .  Theory 

In  Reference  A  (p.  AA).  steady  I  D  flow  (as  in  a  wind  tunnel  lest  section  with 
no  How  distortion)  can  he  represented  by  the  momentum  equation: 


(mass  flow  rate) 

Rr 


(  V  ,  V  ) 

t  Hit  rtlfo 


( 28) 


Drag,  by  definition.  is  a  negative  component  of  VF,  so  swapping  V  and  in 
equation  (28)  will  keep  their  difference  a  positive  value.  The  How  behind  the  cylinder. 


2<i 


for  purposes  in  this  laboratory,  is  two-dimensional.  Rae  &  Pope  develop  the  momentum 
equation  [Ref.  4:p.  214]  into  a  double  integral  and  state  that  "...the  part  of  the  air  that 
passes  over  the  model  suffers  a  loss  of  momentum,  and  this  loss  is  shown  by  and  equal 
to  the  profile  drag..."  Profile  drag,  as  mentioned  earlier,  is  two-dimensional  drag 
composed  of  both  skin  friction  (Df)  and  form  drag  (Dp),  and  is  collectively  denoted  here 
as  D.  Thus,  since  mass  flow  rate  equals  pAV,  where  A  is  the  area  of  the  wake 
perpendicular  to  the  freestream, 


D  -  ffpVjHlV.  -  VJ 


(29) 


In  incompressible  (low  speed)  theory,  pw  equals  pm  (a  constant)  and  the  equation 
becomes 


D  =  ffip.KK  ~  P-O  dA 


(30) 


For  a  unit  span  of  the  cylinder,  as  was  explained  earlier,  S  equals  c  x  1  equals  2r  x  1 , 
and  for  a  unit  slice  of  the  wake  along  the  y-axis,  Da  equals  dy  x  1 .  So  two-dimensional 
drag,  denoted  by  d  may  be  written 


d  -  fipV.yw  ~  P„VW2)  dy  x  1 


(31) 
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Since 


Cn  = 


c  =  d  =  d 

q-(c  *  1}  ’  (±p.Vj)2 r 
2 


(p.yj)r 


and  substituting  equations  (31)  into  (33)  yields 


/  [p.K.y.-p.O  dy 

P.K2r 


The  densities  cancel  out  and  since  the  rest  of  the  denominator  is  also  a  constant,  then 


r  J  v2  V2 


Equation  (9)  from  page  17  may  be  rewritten  as 


la 

-y  ,  while  Vw 


V 


(36) 


Substituting  equation  (36)  into  equation  (35)  yields, 


(37) 


Since  qw  and  qM  are  simply  ((pb)„  -  pj  and  ((Pq).  -  p..),  the  difference  between  these  two 
ratios  are  then  integrated  over  the  width  of  the  wake  to  get 


Similar  results  to  those  in  the  previous  analysis  were  obtained  by: 


1.  Anderson,  [Ref.  2:p.  107,  equation  2.74] 

2.  Rae  and  Pope,  [Ref.  4:p.  215,  equation  4.24] 

3.  Hoemer,  [Ref.  5:p.  1-7] 

4.  Panton,  [Ref.  9:p.  386] 
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rv.  AA2801  CYLINDER  LAB 


A.  PREVIOUS  LABORATORY 

1.  Procedures 

The  former  laboratory  procedures  for  the  AA2801  "Drag  on  a  Cylinder"  may 
be  found  in  Reference  7  and  are  reproduced  in  Appendix  A  for  contrast. 

2.  Limitations  and  Problems 

a.  Limited  Data  Acquisition 

Chief  among  the  limitations  of  the  previous  lab  was  the  lack  of 
numerical  data.  Thus,  all  analysis  was  strictly  drawn  from  graphs  made  by  the  Hewlett- 
Packard  X-Y  Plotters  (see  Figure  12).  Students  would  count  squares  under  the  curves, 
assign  coordinates  to  the  curve,  use  a  planimeter  to  measure  the  area  under  the  curve, 
and  one  student  was  even  known  to  cut  out  the  curves  and  a  small  representative  square 
and  then  weigh  them  on  a  precision  digital  scale!  Also,  not  being  able  to  retrieve  data 
meant  that  this  lab  served  only  one  real  purpose:  Estimating  Drag  on  a  Cylinder. 
Certainly  the  potential  to  expose  students  to  the  functional  use  of  the  wind  tunnel  existed 
but  was  not  exploited  by  this  approach.  Beyond  the  lack  of  data  retrievability,  there 
were  other  problems  with  the  former  lab  procedures. 


29 


Figure  12.  Hewlett  &  Packard  X-Y  Plotters. 


b.  Poor  Relation  Between  Method  Chosen  and  Actual  Procedures 
The  wake  analysis  equation  given  in  the  old  lab. 


Cd  (not  written  as  "cd").  To  begin  with,  the  equation  does  not  agree  with  contemporary 
definitions  of  the  wake  analysis  method  (see  Chapter  HI,  Section  C).  Secondly,  the 
procedure  does  not  support  the  equation  -  either  intuitively  or  precisely.  While  it  was 


pointed  out  that  (V^V,)2  was  simply  qj/q,,  the  V2/V,  portion  was  missing.  Also, 
subtracting  one  curve  from  another  does  not  lead  to  the  solution  of  the  equation  given. 
Finally,  the  P2  -  P,  values  do  not  represent  anything  useful.  P2  was  the  static  source  in 
the  wake  while  P,  was  the  static  port  on  the  cylinder  directed  into  the  freestream  at  6 
equals  0°.  The  difference,  measured  this  way,  was  between  the  static  pressure  through 
the  wake  and  the  stagnation  pressure  at  the  cylinder.  This  was  confusing  and 
bothersome  for  an  otherwise  simple  method  of  drag  estimation. 

c.  Complexity  of  Procedures  and  Sparsity  of  Educational  Value 

To  make  the  plot  described  above,  the  student  had  to  crank  the  traverse 
over  and  back  while  the  technician  had  to  change  leads  into  the  X-Y  plotter.  Graph 
paper  had  to  be  placed  properly  and  the  technician  was  inordinately  relied  upon  to  make 
the  lab  work.  The  student  played  a  minor  role  in  the  lab  other  than  observation  and 
cranking.  Perhaps  one  of  the  more  valuable  parts  of  the  experiment,  the  setting  of  the 
scale  of  the  transducers,  was  left  to  the  technician.  Additionally,  the  student  learned  very 
little  about  data  acquisition,  something  that  would  greatly  benefit  him  or  her  for  future 
use  of  the  wind  tunnel.  The  lab  handout  provided  sparse  information  for  answering  the 
question,  "Why?" 

d.  Inconsistency/Inaccuracy  of  Terminology 

P,  was  never  defined,  nor  was  qj.  HCd"  was  used  for  two-dimensional 
flow  where  "cd"  is  the  convention  used  by  most  texts.  The  length  of  the  cylinder  was 
given  as  28.0  inches  when  it  is  actually  28.4  inches,  which  changes  the  volume  from 
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192.68  in3  to  195.43  in3.  The  pitot  probe  and  the  static  probe  on  the  traverse  were 
referred  to.  as  "two  pitot-static  probes. " 

B.  UPGRADED  LABORATORY 

1.  Data  Acquisition  and  System  Description 
a.  Schematic 

Figure  13  illustrates  the  signal  path  for  the  AA2801  Cylinder  Drag  Lab 
upgrade.  Meanwhile,  the  photograph  in  Figure  14  is  significant:  data  acquisition,  data 
storage,  data  retrieval  and  data  reduction  -  all  in  one  place. 
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Figure  13.  Signal  Path  of  Data  Acquisition  System 
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Figure  14.  486DX/33  MHz  PC  Used  For  Data  Acquisition 
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b.  Signal  Configuration 


(1)  p^,  -  pm.  There  are  two  static  pressure  ports:  a  local  pressure 
port  (pcyl)  located  on  the  surface  of  the  cylinder  approximately  mid-height  and  a  static 
pressure  port  (pw)  located  on  a  traverse  probe  mount.  Both  pressures  feed  into  a  Statham 
pressure  transducer  rated  at  ±  1  psid,  15  Vdc  maximum  input.  The  transducer  is  excited 
by  approximately  11,  volts  and  the  output  signal  is  then  fed  to  an  NPS-built  signal 
conditioner,  and  then  on  to  the  PC. 

(2)  (Po)w  -  Pw  equals  q^  Traverse  mounted  pitot  and  static  probes 
((Po)w  and  pw)  are  co-located  behind  the  cylinder  and  the  pressure  measured  from  these 
are  fed  into  a  Statham  transducer,  Model  P69-1 D-350  rated  at  ±1  psi.  The  signal  path 
to  transducer  and  to  the  PC  is  identical  to  that  of  the  previous  transducer. 

(3)  y.  A  mechanical  traverse  operated  by  a  hand-crank  moves  the 
mount  which  holds  the  pitot  and  static  probes  behind  the  cylinder.  A  small  gear  drives 
a  10-tum,  50K  Q  IRC  potentiometer,  Type  H115-T-5620.  The  potentiometer  is  excited 
with  approximately  11  volts  and  the  output  signal  is  fed  from  the  NPS-built  signal 
conditioner  to  the  PC. 

(4)  0.  The  cylinder  is  mounted  to  a  turntable,  with  rotating  gear 
driving  an  Amphenol  100  Q,  10-tum  potentiometer,  Type  4301  B.  The  potentiometer 
is  excited  by  0.56  volts  while  the  output  signal  is  fed  to  a  Pacific  Amplifier/Signal 
Conditioner,  Model  8255  and  amplified  100  times  before  routing  to  PC.  The  location  of 
the  turntable’s  azimuth  scale  is  such  that  the  cylinder  local  pressure  port  is  pointing 
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directly  into  the  freestream  when  6  equals  0°.  The  turntable  rotates  the  cylinder 
clockwise  from  6  =  -20°  to  8  —  +180°.  Figure  15  shows  the  rotatable  base. 


Figure  15.  Electrically  Rotatable  Base  for  Cylinder 

(5)  (Po)x  -  p„  equals  qx.  The  tunnel  airspeed  indicator  pitot-static 
tube,  located  under  the  ceiling  in  the  forward  part  of  the  test  section,  is  used  to  provide 
(p0)„  and  p„  to  a  Statham  pressure  transducer  rated  at  ±  1  psid,  15  Vdc  maximum  input. 
This  transducer  is  excited  by  approximately  1 1  volts  and  returns  its  output  to  an  NPS- 
built  signal  conditioner  (shown  in  Figure  16),  then  to  the  PC.  The  airspeed  indicator 
serves  as  a  reference  static  pressure  under  the  assumption  that  p«  read  there  remains 
essentially  constant  at  a  fixed  Ap,  as  was  verified  by  calibration. 
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Figure  16.  N.P.S.  Signal  Conditioner 


(6)  CIO  37-pin  MINITERM  Connector  Board.  All  signals  arrive 
here  first  and  are  attached  via  "channel"  points  to  a  standard  37-pin  RS-232  cable  (see 
Figure  13). 

(7)  ComputerBoards  CIO-DAS16/330  A/D  Card.  A  high-speed 
analog  and  digital  data  acquisition  card  was  placed  in  the  PC  with  a  base  address  of 
300H  (300  HEX).  A  special  library  in  the  software  package  allowed  the  computer  to  use 
the  A/D  card  (see  Appendix  B). 


(8)  IBM  486DX/33  MHz  PC.  A  standard  486  IBM-clone  PC  was  the 
terminal  for  the  data  acquisition  system.  This  PC  had  DOS  version  S.O  and  Windows 
version  3.1  loaded. 

c.  Software 

LabTech  Notebook  for  Windows,  version  7.1.1,  was  the  operating 
software  used  to  drive  the  data  acquisition  system.  It  is  a  process  and  control  type 
application  which  can  be  used  to  both  collect  data  and  to  control  systems. 

LabTech  Notebook  has  two  basic  operating  modes:  BUILD-TIME  and 
RUN-TIME.  As  their  titles  imply,  BUILD-TIME  is  used  to  construct  set-ups  which  are 
then  executed  in  RUN-TIME.  These  set-ups  contain  customized  "ICONblocks"  which 
connect  signals  and  create  displays  and  data  files. 

2.  "ZERO"  and  "SPAN" 

These  two  LabTech  Notebook  set-ups  were  created  to  handle  the  job  of 
scaling  the  signals  coming  into  the  program.  It  was  felt  that  the  signals  must  be  correct, 
not  just  leaving  the  signal  conditioners,  but  more  importantly,  entering  the  PC.  Noise 
and  other  anomalies  affect  the  signal  at  all  points  in  its  path.  ZERO  and  SPAN  are  fully 
explained  in  Appendix  B. 

3.  Procedures 

a.  Run  #1  -  Pressure  Distribution  Method 
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(1)  Approach.  While  maintaining  the  reference  static  pressure  pw 
outside  the  wake  and  behind  the  cylinder  (which  is  thus  pM),  the  cylinder  is  rotated  at  a 
constant  rate  (and  so  is  pcyl)  from  0°  to  +180°.  The  turntable  is  started  at  0  =  -20°  to 
ensure  any  acceleration  is  essentially  zero  by  the  time  it  passes  0  =  0°  and  the  maximum 
pcy,  is  actually  recorded,  although  this  is  not  required. 

(2)  Measured  sources.  The  signals  being  measured  were  pcy,  -  pw  v. 
0,  where  pw  is  the  reference  pressure  and  is  assumed  to  be  the  freest  ream  static  pressure 
(Po.)  of  the  test  section  and  where  0  is  the  position  of  pcy„  the  local  static  pressure  on  the 
cylinder,  as  it  varies  from  the  stagnation  pressure,  through  p,,  past  the  point  of 
maximum  velocity,  past  the  separation  point  of  the  boundary  layer  and  finally  to  0  equals 
180°. 

(3)  Purpose.  These  measurements  will  be  used  by  the  software  to 
find  and  plot  Cp  v.  0. 

(4)  Signal  scaling  criteria.  Before  running  the  tunnel,  at  zero  tunnel 
Ap,  since  q,  pcy„  p„,  pw  and  Cp  all  equal  0,  the  transducer  output  is  conditioned  to  also 
read  0.00  volts.  Next,  the  tunnel  is  run  at  15  cm.  H20,  with  0  equal  to  0°  and  y  equal 
to  0.0  inches.  But  at  0  equals  0°,  pcyl  equals  (p0)«..  Since  it  is  known  that 
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(Po)o.  equals  q„  +  p„,  the  following  substitution  can  be  made  for  Cp: 


at  0  =  0°,  C 

P 


(<?.  *  PJ  ~  Pm 
<7- 


(40) 


Knowing  this,  the  output  signal  can  be  conditioned  for  the  transducer  to  read  0.01  volts. 
The  software  used  on  the  PC  is  then  simply  configured  to  multiply  the  signal  coming  in 
by  100  so  that  0.01  volts  is  displayed  as  1.0  when  Cp  actually  equals  1.0.  Figure  17 
illustrates  the  LabTech  Notebook  set-up  "PRE  DISl."  Figure  18  shows  the  actual 


Figure  17.  LabTech  Notebook  Set-Up  "PRE_DIS1” 
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results  of  executing  PRE  DISl  in  RUN-TIME. 

(5)  Expected  results.  From  inviscid,  incompressible  theory,  the 
pressure  distribution  method  will  generate  a  cosine-type  curve  starting  out  with  Cp  equal 
to  1 .0  at  0  equal  to  0°  and  decreasing  through  Cp  equal  to  0  to  a  minimum  of  Cp  equal 
to  -3.0  or  separation,  whichever  occurs  first  (see  Figure  18). 


Figure  18.  Actual  LabTech  Notebook  Readout  of  Cp  vs.  0  from 
Run  n\ 


The  point  of  separation  will  also  define  the  point  of  wake  initiation.  It  is  anticipated  that 
the  wake  behind  a  stationary  cylinder  in  moving  fluid,  where  the  Reynolds  numbers  are 
of  the  order  2.7  x  10s,  will  not  change  appreciably  in  width  (in  the  y-direction)  from 
where  it  begins  [Ref.  2:p.  238,  Fig.  3.40 d\.  As  the  wake  progresses  in  the  x-direction 
away  from  the  cylinder,  it  seeks  to  recover  back  to  the  freestream  condition.  Thus,  the 
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width  of  the  wake,  for  this  Reynolds  number,  may  be  approximated  by  entering  the  value 
for  azimuth  at  the  point  of  separation  into  the  following  equation: 


maximum  expected  wake  width  (y^,)  =  2  (  d  |  sin  0  |  ) 
where  d  is  the  diameter  of  the  cylinder  in  inches.  This  value  was  used  to  predict  the 
beginning  of  the  wake  during  Run  #2. 

(6)  Using  the  data  obtained  during  data  acquisition  to  estimate  the 
drag  coefficient.  Using  the  pressure  distribution  method  for  blunt  bodies,  form  drag  is 
the  significant  contributor  to  the  resultant  force  on  a  non-rotating  cylinder  in  airflow 
[Ref.  3:pp.  59,  60],  [Ref.  2:p.  60,  209]: 


C  m  f*  cos6  60,  where  C„  =  ^ 

'  Jo  9.  '  9„ 


(42) 


The  data  obtain  during  Run  #1  is  a  column  of  raw  Cp  values,  a  column  of  corresponding 
azimuth  (0)  values,  and  a  column  of  theoretical  Cp  values,  the  latter  of  which  are 
obtained  from  the  equation  for  theoretical  inviscid,  incompressible  surface  pressure 
coefficients  over  a  cylinder  [Ref.  2:p.  199]: 

C  =  1  -  4sin20  (43) 

p 


Having  these  data,  any  program  or  technique  of  integrating  the 
first  two  columns  over  the  interval  of  0  to  t  will  supply  the  value  cdp.  Tabular  data  were 
not  available  in  previous  labs  and  students  relied  on  physically  measuring  the  area  under 
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the  curve  from  0°  to  +180°.  Now,  using  LabTech  Notebook  and  automatic  data 
acquisition,  data  are  stored  in  an  ASCII  file  and  may  be  imported  into  a  program  such 
as  MATLAB  for  analysis. 

To  compare  the  obtained  from  the  experiment  to  empirical 
data,  the  Reynolds  number  is  calculated  from 


Re  * 


pVJ 


(44) 


and  the  graph  of  Re  vs.  CD  shown  in  Reference  7,  (p.  36)  or  Figure  24  is  used.  Care 
must  be  taken  to  correct  Re  for  the  turbulence  factor,  TF,  to  produce  the  effective  Re: 

^eEFF  ~  ^enncorrccuJ^ 

where 

C  =  +  (46) 


b.  Run  #2  -  Wake  Analysis  Method 

(1)  Approach.  While  the  cylinder  is  fixed  in  some  (any)  position  6, 
the  wake  probes  are  moved  along  a  y-axis  traverse  from  0.0  inches  to  +20.0  inches. 
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(2)  Measured  sources.  Three  values  are  measured  during  this  run: 
qw,  q,,  and  traverse  location,  y.  By  definition,  (p0)w  -  pw  equals  qw  and  (p0)o,  -  p„  equals 
q„.  It  has  been  shown  experimentally  that  when  the  probes  which  measure  (p0)w  -  pw  are 
located  outside  the  wake,  qw  equals  q».  Since  qw/qw  is  the  ratio  sought,  the  wake 
analysis  (momentum)  method  is  really  used  to  determine  variations  in  velocity  through 
the  wake. 

(3)  Purpose.  These  values  will  be  used  by  the  software  to  plot  the 
following  difference  in  ratios: 


which  is 


\|  (Po)-  -  P. 


top)*  p"  v 
(Po>-  -  P- 


(47) 


(48) 


The  difference  between  these  two  ratios  forms  the  integrand  of  the  momentum  method 
for  calculating  the  drag  coefficient,  as  was  explained  earlier. 

(4)  Signal  scaling  criteria.  Before  running  the  tunnel,  at  zero  tunnel 
delta  P,  all  pressures  are  equal,  all  q’s  are  zero  and  all  transducer  outputs  are  scaled  to 
read  0.0  volts.  Next,  the  tunnel  is  run  at  15  cm  H20,  with  y  equal  to  0  inches.  Here, 
outside  the  wake,  qw  equals  q„  (both  A  0.0).  Since  q^/q*  will  be  plotted,  using  the 
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Commutative  Law,  both  outputs  are  scaled  to  the  same  value,  in  this  case  1 .0.  Actually, 
any  common  value  would  have  worked  since  the  ratio  when  qw  equals  q*  is  always  unity. 
Figure  19  shows  the  set-up  WAKEANl/2  from  LabTech  Notebook. 


Figure  19.  LabTech  Notebook  Set-Up  "WAKE_ANI/2" 


(5)  Expected  results.  From  the  value  of  maximum  expected  wake 
width  (y’w,kc)  estimated  during  Run  H\  (equation  (43)),  no  appreciable  change  in  pressure 
should  be  seen  until  reaching  the  interval 


y»ak*  inches)  -  10  ± 


/ 

y^akt 


(49) 
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Logically,  q*  should  remain  constant  during  the  tunnel  run.  However,  qw,  the  local 
dynamic  pressure  in  the  wake,  should  change  to  reflect  the  condition  of  the  wake  as  the 
traverse  moves  from  0.0  to  20  inches.  Also,  pw  should  change  less  than  (po)w.  Indeed, 
many  "profile  drag"  rakes  used  to  measure  wake  drag  from  attached  flow  over  airfoils, 
often  have  only  a  single  static  probe  for  dozens  of  total  pressure  probes  [Ref.  4:p.  1 12]. 
The  static  pressure  is  treated  to  essentially  remain  constant  through  the  wake.  Actually, 
due  to  tunnel  wall  effects,  the  static  pressure  may  change  slightly  [Ref.  l:p.  35]. 
Nevertheless,  it  is  expected  that  as  the  wake  is  entered,  qw  should  get  smaller  because 
(p0)w  decreases  at  a  faster  rate  than  pw. 


(6)  Using  the  data  obtained  during  data  acquisition  to  estimate  the 
drag  coefficient.  The  software  set-up  "WAKE  AN1 "  in  LabTech  Notebook  is 
configured  to  receive  the  three  inputs:  q»,  q„  and  y.  Through  a  string  of  "calculation 
blocks,"  WAKE  AN1  forms  the  plot  mentioned  earlier  of 


(Pp)w  P  W 

^  (Po).  -  P. 


(Pp)w  -  Pw  v 
<Po).  ■  P- 


which  is 


(50) 


(51) 
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The  ordinate  of  the  plot  is  the  integrand  from 


c  =  1  f20  (  (Pq)w~Pw  _  (P()W~PW  .  (52) 

*  r  Jo  ^  <pjm-pm  (Po)m~Pm 

t 

This  equation  is  integrated  from  0.0  to  20.0  inches  and  the  result  is  an  approximation  of 
the  drag  coefficient.  Although  the  setup  WAKE_AN1  creates  a  display  of  the  ratios  of 
dynamic  pressure  versus  traverse  distance  y,  three  columns  of  data  are  also  generated  for 
the  student  to  analyze,  using  a  program  such  as  MATLAB.  They  are  simply  the  three 
inputs  qw,  q.  and  y.  It  is  a  simple  matter  to  write  a  code  to  use  these  inputs  to  solve 
for  cd.  See  Appendix  D  for  an  example. 

c.  Run  #3  -  Wake  Analysis  Method  with  Tripped  Boundary-Layer 

The  only  difference  between  Run  it 2  and  Run  #3  is  the  addition  of  two 
"cellophane"  adhesive  tape  strips,  Vi"  x  20"  x  1  mil.  These  strips  are  affixed  to  the 
cylinder,  approximately  1.25"  on  either  side  of  the  0  equals  0°  position,  parallel  to  the 
z-axis  (vertical  axis). 

The  function  of  these  two  strips  are  to  act  like  tiny  "vortex  generators," 
which  are  familiar  to  all  Naval  aviators.  The  strips  effectively  "trip"  the  attached 
boundary  layer  causing  otherwise  laminar  flow  to  become  turbulent.  The  consequence 
of  this  act  is  dramatic.  By  inducing  an  early  transition  from  laminar  to  turbulent  flow, 
the  engineer  can  delay  actual  flow  separation,  thus  delay  the  formation  of  the  wake.  This 
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in  turn  decreases  the  drag  produced  by  an  object  in  an  airflow,  a  favorable  situation  for 
blunt  bodies  [Ref.  2:p.  644]. 

There  is  an  intuitive  explanation  for  why  flow  separation  is  delayed  by 
creating  turbulence  where  there  once  was  laminar  flow:  In  the  very  small  boundary 
layer,  laminar  flow  can  be  thought  of  as  flow  in  1  dimension  only  -  parallel  to  the 
surface.  There  are  no  significant  forces  at  work  to  move  the  flow  either  closer  toward 
the  surface  or  farther  away.  Eventually,  however,  the  displacement  of  the  flow  from  its 
original  freestream  position  becomes  so  great  that  the  laminar  boundary  layer  flow  leaves 
the  surface  of  the  airfoil  and  seeks  to  return  to  its  former  place.  If,  however,  this 
laminar  flow  is  artificially  brought  into  contact  with  the  surface  of  the  airfoil,  which  is 
by  no  means  smooth,  the  air  molecules  will  be  caused  to  trip  and  tumble.  Specifically, 
in  some  random  direction  in  relation  to  the  surface  of  the  airfoil,  but  in  general  in  a 
continued  downstream  direction.  As  this  air  trips  and  tumbles  (turbulence)  and  moves 
along  the  surface  of  the  airfoil,  it  acts  like  a  cartoon  snowball  rolling  down  a  winter 
hillside.  As  it  grows,  it  "reaches"  farther  above  the  surface  of  the  airfoil  and  "gathers" 
in  higher  energy  air.  This  higher  energy  causes  the  air  near  the  surface  to  become  much 
faster  than  it  was  when  it  was  laminar.  This  new  speed  in  turn  causes  the  flow  to  "push” 
ahead  farther  along  the  surface  than  it  naturally  wants  to  go.  Eventually  the  forces  at 
work  trying  to  restore  all  flows  to  the  original  freestream  condition  overcome  the  contact 
between  the  now  turbulent  flow  and  the  surface  and  the  airflow  separates  and  is  returned 
to  the  freestream.  Now  the  point  of  separation  is  much  farther  along  and  thus  the  wake 
width  and  associated  drag  is  significantly  reduced. 
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Anderson,  in  Reference  2  (p.  643),  states: 


Because  of  the  agitated  motion  in  a  turbulent  flow,  the  higher-energy  fluid 
elements  from  the  outer  regions  of  the  flow  are  pumped  close  to  the  surface. 
Hence,  the  average  flow  velocity  near  a  solid  surface  is  larger  for  a  turbulent 
flow  in  comparison  with  laminar  flow... because  the  energy  of  the  fluid 
elements  close  to  the  surface  is  larger  in  a  turbulent  flow,  a  turbulent  flow 
does  not  separate  from  the  surface  as  readily  as  a  laminar  flow.  If  flow  over 
a  body  is  turbulent,  it  is  less  likely  to  separate  from  the  body  surface,  and  if 
flow  separation  does  occur,  the  separated  region  will  be  smaller.  As  a  result, 
the  pressure  drag  due  to  flow  separation,  Dp  [form  drag],  will  be  smaller... 


Though  not  demonstrated  in  this  lab,  Hoemer  points  out  in  his  classic 
Fluid  Dynamic  Drag  (p.  3-27)  yet  another  way  to  reduce  separation  around  a  stationary 
cylinder:  use  suction  on  the  cylinder  via  a  porous  surface.  This  prevents  or  delays  the 
start  of  an  alternating  vortex  sheet. 

The  helicopter  industry  is  already  producing  aircraft  which  use  a 
delayed-separation  principle  called  the  "Coanda"-effect.  The  application  is  NOTAR”, 
which  stands  for  "No-Tail-Rotor"  [Ref.  10:p.  If].  Here  blowing,  vice  suction,  is  used 
to  manipulate  the  point  of  separation  around  a  cylindrical  tail  boom  to  provide  an 
asymmetric  surface  pressure  distribution,  which  then  generates  desired  lift. 
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V.  RESULTS 


A.  ACTUAL  CLASS  DATA  (USING  DRAG1.M) 

Table  5  is  a  summary  of  the  results  obtained  by  using  the  program  dragl  .m  created 
in  Matlab  code,  which  is  located  and  explained  in  Appendix  D.  Users  of  dragl  .m  may 
access  the  explanation  by  typing 

>  >help  dragl 

Dragl  m  assumes  the  following:  NPS  Low  Speed  Wind  Tunnel  and  associated  TF 
and  calibration  factor;  tunnel  operated  at  15  cm  H20;  equations  contained  in  this  thesis 
are  used;  data  have  been  edited  to  remove  "idle  data.  ” 

Table  5  results  are  from  actual  classroom  data  obtained  from  the  three  LabTech 
Notebook  set-ups  described  in  Chapter  IV.  Two  AA2801  labs  of  approximately  seven 
students  each  were  conducted  using  the  upgraded  laboratory  handout  contained  in 
Appendix  B.  Since  this  lab  was  conducted  on  9  November  1993,  Appendix  B  has  been 
revised  to  include  the  lab  times  required,  Table  1  -  Drag  Classification,  and  illustrations 
under  the  ZERO/SPAN  paragraph. 
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TABLE  5.  DRAGl.m  results 


»  dragl 

§  %  (from  data  taken  during  9  Nov  1993  %  1 

|  AA2801  Lab) 

1  % 

%  This  data  has  been  reduced  using 

%  DRAGl.m.  It  has  been  corrected 

%  for  blockage  and  is  explained  in 

%  APPENDIX  D 

% 

cdpdu 

=  0.71 

cdpdl 

=  0.68 

cdpd2 

=  0.67 

wakebeg 

=  7.3  in 

wakeend 

=  12.7  in 

sepaz 

«  64.1° 

pftl 

=  0.68 

Reeff 

*  2.89e+05 

cdwal 

=  0.99 

cdwa2 

=  0.18 

percredu 

=  81.7  % 

B.  PRESSURE  DISTRIBUTION  METHOD  (RUN  #1) 

1.  Data 

Data  were  obtained,  as  explained  in  Appendix  B,  and  filed  under 
PREDISl  .m.  Actual  results  from  the  classroom  are  recorded  in  Appendix  C.  It  should 
be  noted  that,  in  an  effort  to  save  space,  this  DOS  Text  File  has  been  saved  in  small 
print  and  placed  into  multiple  columns. 

The  set-up  PRE  DISl  automatically  begins  collecting  data  when  the  rotator 
passes  -17°  and  stops  when  it  passes  +179°.  This  worked  because  the  rotator  itself 
turned  at  a  constant  rate  and  was  controllable  electrically.  The  data  were  originally 
saved  in  two  columns:  0  and  Cp(raw).  All  data  prior  to  0°  were  manually  edited  out. 
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Due  to  the  assumed  constant  rate  of  the  rotator  combined  with  the  constant 
sampling  rate  of  10  Hz,  the  trapezoidal  rule  was  applied  to  the  data  using  [Ref.  8:p. 
193]: 

(53. 

An  uncorrected  (for  blockage  factors)  c^  of  0.71  was  obtained  using  a 
modified  trapezoidal  rule: 

(54> 

Here,  equispaced  x-values  were  not  assumed. 

The  blockage  factors  described  earlier  (equations  (3)-(5))  were  applied  to  the 
results  from  both  methods,  yielding  a  c^  of  0.68  using  equation  (55)  and  0.67  using 
equation  (56).  Figure  20  shows  how  the  MATLAB  program  Dragl  .m  plots  Cp  vs.  0,  and 
places  an  asterisk  where  separation  is  believed  to  occur.  A  routine  which  assumed  Cp 
increases  abruptly  during  classical  cylinder  flow  separation  predicted  the  separation 
azimuth  to  be  64. 1  °.  This  value  for  theta  was  then  used  with  equations  (43)  and  (51)  to 
estimate  the  wake  width,  which  ranged  from  7.3  to  12.7  inches. 

Next,  Matlab-defined  functions  "polyfit"  and  "polyval”  were  used  to  generate 
a  fifth-order  polynomial  approximation  to  the  actual  data.  In  this  case,  a  fifth-order 
polynomial  was  chosen  as  it  best  fit  the  plot  given.  Results  are  shown  in  Figure  21 . 
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Figure  20  The  Separation  Point  Estimated  by  SEPAZ 


Finally,  the  effective  Reynolds  number  of  2.89  x  10'  was  calculated  using 
equations  (7),  (8),  (13)  and  (46). 

2.  Analysis 

The  corrected  c.,  s  plot  close  to  but  not  on  the  empirical  data  shown  in 
Figure  24.  One  reason  why  0.67  is  seemingly  low  is  that  the  pressure  distribution 
method  does  not  measure  c,t;  it  only  measures  the  component  due  to  pressure  (form 
drag).  In  fact.  Panton,  in  Incompressible  Flow  [Ref.  9:p.  386|.  was  responsible 
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C.  WAKE  ANALYSIS  METHOD  (RUNS  #2.3) 

I.  Data 

Essentially  the  same  methods  just  explained  were  used  to  find  c.,  related  to 
the  wake  analysis  from  Chapter  V.  Section  B.  Equation  (56)  was  used  to  find  ct,  for 
Run  #2  (0.99)  and  Run  #3  (0.18).  The  wake  analysis  method  uses  the  WAKE  AN  1/2 
set-ups  and  from  data  reduction  using  DRAGI.m,  the  following  plots  were  generated 


(Figures  22  and  23): 
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of  the  data  files.  These  idle  data  were  easy  to  detect.  All  data  prior  to  the  y-value  (0) 
increases  and  all  data  after  it  stops  increasing  were  edited  out. 

Figure  24  shows  the  results  of  all  three  runs  plotted  against  emperical  data. 


Figure  24.  Results  Plotted  Against  Emperical  Data 
2.  Analysis 

Again,  the  corrected  cd  from  Run  #2  was  plotted  against  the  Reefr  and 
compared  to  the  empirical  data.  It  compared  very  favorably,  tending  to  validate  the  TF 
given.  Also,  from  Figure  22,  the  wake  analysis  width  prediction  of  DRAGl.m  seemed 
reasonable. 

The  81 .7%  reduction  in  cd  from  Run  #3  is  dramatic.  If  this  number  is  high, 
then  the  TF  might  come  into  question  again.  "The  higher  the  freestream  turbulence,  the 
more  readily  transition  takes  place  [Ref  2 :p.  389] . "  Golf  balls  use  dimples  to  accomplish 
similar  levels  of  dramatic  drag  reduction. 
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There  is  one  concern  that  arises  out  of  Figures  22  and  23.  It  would  appear 
that  given  the  use  of  equation  (53)  it  should  be  quite  impossible  to  ever  have  a  negative 
solution,  as  was  seen  during  both  runs.  There  are  two  possible  explanations  for  this 
phenomenon:  (1)  As  shown  in  Figure  25  and  as  discussed  on  page  35  of  Reference  1 ,  it 
is  theoretically  possible  to  have  qw  become  larger  than  q„.  Were  this  to  happen,  the  plot 
of  equation  (53)  would  indeed  be  negative  for  that  interval.  (2)  A  more  likely  reason, 
though,  is  that  the  signal  from  the  qw  or  q^  transducer  is  either  poorly  scaled  or  actually 
drifts  with  either  time  or  with  displacement.  This  could  be  easily  verified  by  repeating 
a  traverse  of  the  test  section  in  both  directions,  with  the  cylinder  removed,  and 
comparing  qw  to  q„.  Unfortunately,  time  did  not  allow  for  this  verification  and  a 
solution  to  the  negative  values  seen  in  Figures  22  and  23  are  left  for  future  exploration. 
For  this  thesis,  the  data  were  neither  adjusted  up  nor  down  (from  the  zero  axis),  nor  was 
the  plot  in  Figure  22  pivoted  by  anchoring  the  left  end  and  rotating  the  right  end  of  the 
plot  up  to  the  zero  axis.  These  corrective  methods  might  result  in  a  better 
approximation,  but  it  remains  to  be  seen  whether  data  acquisition  or  tunnel  velocity 
profiles  are  responsible  for  the  negative  values. 
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Figure  25  Graph  of  flic  Equation  y-sqrt(x)  x 

VI.  CONCLUSIONS 

A.  CPCRADF  LIMITATIONS 
I.  Wake  Analysis  MefhoH 

Rae  and  Pope  point  out  in  Reference  4  (pp  214  21 7),  that  the  wake  analysis 
method  lias  several  "gray  areas:" 

1  How  do  you  measure  spanwise  variation  of  the  profile  (hag 
behind  an  airfoil?  (p.  213) 

2  Are  vortices  present  at  higher  Re?  |p  21 4| 


+  G  ) 

(see  Ref.  9:p.  385) 

.  ■  i  ....  - . —  — ...  ».  .  . . 

Tills  plot  shows  then,  (see  Ref.  1:p.  35) 

what  would  happen  if 
q(wake)  =->  q(lnflnlty),  tliat  is, 
if  the  ratio  (qw/qi)  was  ever 
>  1.0  (shown  here  yielding  a  ^ 
negative  result  on  the  interval  / 

1.0  x  2.0)...  / 


SR 


3.  Is  true  two-dimensional  flow  a  function  of  Re  or  does  it  really 
exist?  [p.  214] 

The  authors  state  twice: 

It  is  well  known  that  the  wake  survey  is  not  valid... where  separation  is 
present... cannot  measure  momentum  loss  caused  by  separation  or  fluid 
rotation,  [p.  214,217] 

The  limitations,  then,  are  that  it  is  yet  unknown  exactly  how  valid  the 
approximation  of  cd  is  from  the  wake  analysis  method.  Certainly  flow  separation  is 
present.  While  the  value  obtained,  cd  =  0.99,  fits  that  empirical  data,  enough  questions 
are  raised  by  Rae  and  Pope  that  no  strong  conclusions  may  be  drawn  at  this  time. 

2.  Hardware 

Current  hardware  returns  signals  on  the  order  of  l(f3  volts.  Subsequent 
amplification  of  these  tiny  signals  has  the  effect  of  amplifying  noise  or  simply  limiting 
the  usefulness  of  the  signal.  The  Analog-to-Digital  data  acquisition  card  used  in  this  lab 
has  the  ability  to  resolve  signals  down  to  409S  parts.  In  some  cases,  the  only  signals 
possible  were  confined  to  a  range  of  0.000  to  0.009  volts  (see  Table  3,  Appendix  B). 
Meanwhile,  the  smallest  range  LabTech  Notebook  can  program  down  to  read  is  ±  0.625 
volts.  This  means  that 

*0.625  -  (  -  0.625)  vote  „  ±  aooo33  vote  mobaim  (55) 
(4096  -  1  )  counts  count 

Dividing  an  actual  useful  range  of  0.009  volts  by  this  resolution  yields  an  approximate 
effective  resolution  down  to  just  28  signals  discriminating  between  0.000  and  0.009  volts. 
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Put  another  way,  while  the  A/D  card  can  resolve  down  to  4095  parts,  only  28  of  these 
parts  are  being  used,  less  than  1  %  of  the  card’s  capability. 

3.  LabTech  Notebook 

This  software  application  requires  considerable  time  to  learn.  Having  a 
qualified  user  available  helps  greatly.  The  lengthy  tutorial,  once  mastered,  provides 
an  excellent  foundation  for  using  this  software  package  to  do  data  acquisition. 

Beyond  the  time  it  took  to  leant,  it  has  three  other  limitations: 

1.  apparent  limited  sampling  rate  of  1000  Hz 

2.  very  limited  graphics  capability 

3.  high  memory  usage  in  the  Windows  environment. 

B.  RECOMMENDATIONS 

1.  Study  the  Wake  Analysis  Method  Further 

A  survey  of  the  static  pressure  in  the  wake  behind  the  cylinder  should  be 
taken  and  compared  to  the  freestream  static  pressure.  If  the  static  pressures  vary 
significantly,  the  cylinder  should  be  replaced  and  the  tunnel  velocity  profile  measured 
along  the  traverse,  paying  close  attention  to  the  ends.  The  velocity  along  the  traverse 
should  be  checked  for  speeds  which  are  greater  than  freestream.  Rae  and  Pope  warn 
that  q  must  be  taken  as  far  behind  the  cylinder  so  that  pw  equals  p*  [Ref.  4:p.  216],  If 
the  wake  velocity  remains  less  than  the  freestream  velocity,  check  for  signal  noise  or 
other  faulty  modes  of  the  data  acquisition  system. 
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2.  Match  Hardware  to  Software  Capability 

The  12-bit  resolution  of  the  A/D  card  as  just  explained,  combined  with  as 
large  an  output  as  possible,  would  yield  a  higher  resolution  than  currently  available. 
This  means  that  the  signals  into  the  PC  should  have  as  wide  a  range  as  possible,  not  the 
current  9  millivolt  range  output  currently  seen.  This  output  is  both  a  function  of  the 
signal  transducer  and  potentiometer  and  signal  conditioner.  If  the  hardware  used  to 
acquire  signal  sources  must  remain,  then  the  signal  conditioner  must  be  reconfigured  to 
amplify  the  output.  If  the  sources  can  be  measured  using  devices  which  discriminate 
well  below  .033  millivolts  level,  then  the  signal  conditioners  need  not  be  changed. 

3.  Investigate  Better  Data  Acquisition  Software 

Investigate  non-Windows  applications  for  better  response  to  memory 
limitations  of  PC.  Also,  recommend  obtaining,  for  long-time  use,  a  data  acquisition 
program  which  has  a  broader  graphics  package  than  LabTech  Notebook  currently 
employs. 

4.  Validate  TF 

Since  no  attempt  was  made  to  validate  the  turbulence  factor,  the  TF  should 
be  recalibrated  using  procedures  discussed  in  Reference  1.  TF  is  important  and 
contributes  approximately  equally  with  the  blockage  factors  in  the  calculation  of  the 
effective  Reynolds  number  (equation  (8)). 

5.  Upgrade  the  NPS  Low  Speed  Wind  Tunnel  deck 

The  next  logical  step  in  the  wind  tunnel  upgrade  is  to  match  the  data 
acquisition  capability  with  the  actual  layout  of  the  tunnel  test  facilities.  Layers  of  old 
thesis  projects,  built  upon  each  other  over  the  years,  can  effectively  be  replaced  with  this 
data  acquisition  system.  Former  procedures  and  techniques  may  be  incoiporated  into  the 
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set-ups  in  LabTech  Notebook  on  an  as-needed  basis.  Furthermore,  the  current  collection 
of  signal  conditioners,  signal  amplifiers,  multiplexer  and  wiring  inhibits  easy  access  to 
the  test  section.  The  PC  used  in  this  thesis  must  be  located  a  full  twenty  feet  away  from 
the  test  section  making  lab  observation  and  display  observation  all  but  exclusive  events. 
A  permanent  wood  lab  table  extending  from  the  test  section  along  the  diffuser,  would  be 
a  practical  way  to  provide  this  easy  access  and  give  a  place  for  lab  equipment,  in 
component  form,  to  be  set  up.  Furthermore,  a  system  of  annually  justifying  the 
continued  presence  of  former  thesis  projects  is  required  to  prevent  the  problem  of 
layering.  Copies  of  thesis  reports  should  be  readily  available  near  each  project  left  intact 
after  the  author,  student  or  professor,  transfers  from  the  school. 
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APPENDIX  A:  PREVIOUS  LAB  HANDOUT 


Drag  on  a  Cylinder 


Objectives 

The  drag  on  a  circular  cylinder  which  spans  the  height  of  the  wind  tunnel  test  section 
will  be  deteimined  using  wake  analysis  and  pressure  distribution  methods,  demonstrating 
basic  wind  tunnel  methodology,  including  operating  procedures,  Reynolds  number 
calculation  and  application  of  wind  tunnel  correction  factors. 

Apparatus 

The  NPS  low  speed  wind  tunnel  is  described  in  detail  in  the  NPS  Laboratory  Manual  for 
Low  Speed  Wind  Tunnel  Testing.  A  circular  cylinder  (length  =  28.0  in.,  diameter  = 
2.96  in.,  volume  =  192.68  in3)  is  placed  in  the  wind  tunnel  such  that  it  spans  the  test 
section  height.  A  single  static  port  is  built  into  the  mid-portion  of  the  cylinder.  The 
cylinder  itself,  is  attached  to  an  electrically  operated  rotatable  base.  In  the  aft  portion 
of  the  test  section  are  placed  two  pitot-static  probes,  aligned  such  that  one  provides  total 
pressure  and  the  other  static  pressure,  in  the  same  longitudinal  plane.  These  probes  are 
attached  to  a  traverse  mechanism  which  enables  lateral  positioning  (across  the  width  of 
the  test  section)  to  be  varied.  Pitot-static  pressures  are  transduced  by  strain  gauges  and 
output  voltages  fed  to  x-y  plotters  (y-information).  Additionally,  sensors  on  both  the 
rotatable  base  and  traverse  mechanism  provide  position  signals  to  the  plotters  (x- 
information).  Thus,  wind  tunnel  output  consists  of  five  electrical  signals: 

Pwake  ~  ?! 

Qwake  9 2 

P cylinder 

Cylinder  position  (6) 

Probe  location 

The  tunnel  will  be  run  at  a  pressure  differential  of  15  cm  H20  (approximately  1 15  kts.) 
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Wake  Analysis  Methodology 

Applying  the  Reynold's  [sicj  Transport  Theorem  to  Newton’s  Second  Law, 

£  F  =  —  ( momentum ) 

6t 


and  writing  the  resultant  equation  for  a  control  volume  (see  Zucker:  Fundamentals  of  Gas 
Dynamics,  p.  65-69)  results  in: 


C  - 

d  d"vx 


(4.1) 


Pressure  Distribution  Methodology 

The  forces  on  a  cylinder  may  be  modeled  as  shown  in  Figure  2 1 . 


Below  the  critical  Reynold’s  [sic]  number,  the  form  drag  is  the  dominant  component  of 
the  resultant  force  on  the  cylinder.  (Typically,  total  parasite  drag  is  composed  of 
approximately  80%  form  drag  and  20%  skin  friction  drag.)  For  blunt  bodies  (see  Bertin 
and  Smith:  Aerodynamics  for  Engineers,  p.  88-91),  form  drag  can  be  approximated  by 
the  expression: 
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Cd~f*cpcosQ{>6 


COS050 


(4.2) 


X-Y  Plotter  Calibration 

X-Y  plotters  can  be  calibrated  to  provide  the  following  ratios: 

-  =  (YJV,)1  vs.  y 


*1 


Vi.  y 


P  -  P 
*6  r\ 


vs.  0 


From  the  curves  of  these  ratios,  the  integrals  of  equations  (4.1)  and  (4.2)  may  be 
evaluated.  Calibration  of  the  plotters  is  accomplished  as  follows: 

.  (VJVJ2  vs.  y: 

1 .  The  plotter  pen  is  positioned  at  the  zero  point  prior  to  starting  the  tunnel  (V2  =  0). 

2.  After  the  tunnel  is  at  full  speed,  with  probes  located  away  from  the  wake  (V2  = 
V,),  the  plotter  pen  is  positioned  at  the  (V2/V,)2  =  1.0,  y=0  point. 

3.  The  probes  are  then  traversed  across  the  test  section  (20  in.)  and  appropriate  ratios 
are  plotted. 
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vs.  y : 


P2  ~Pi 

•  - 

1 .  After  the  tunnel  is  at  full  speed,  with  probes  located  away  from  the  wake  (P2  = 
P,),  the  plotter  pen  is  positioned  at  the  Cp  =  0.0,  y=0  point. 

2.  The  probes  are  traversed  across  the  test  section  and  the  appropriate  ratios  versus 
y  are  plotted  simultaneously. 


*1 

1 .  The  plotter  pen  is  positioned  at  the  Cp  =  0  point  prior  to  starting  the  tunnel. 

2.  With  the  tunnel  running,  the  cylinder  is  rotated  such  that  0  =  0°,  which  puts  the 
static  port  at  the  stagnation  point  (P2  =  P,  +  q,,  and  the  ratio  equals  unity).  The 
plotter  pen  is  positioned  at  the  Cp  =  1.0,  6  =  0°  point. 

3.  The  cylinder  is  rotated  to  the  8  =  -20°  position  and  plotting  begins  (from  8  =  -20° 
to  180°). 

Laboratory  Procedures 

1 .  Determine  atmospheric  pressure: 

patm  ~  (baro.  reading) + (temp  corr.)+(lat.  corr.) 

(Latitude  correction  =  -.0245) 

2.  Determine  temperature. 

3.  Make  required  wind  tunnel  runs. 

4.  Trip  the  boundary  layer  and  make  an  additional  run  in  order  to  determine  the 
reduction  in  wake  size  afforded  by  a  turbulent  boundary  layer.  (Wake  analysis 
methodology  only.) 

Conversions  and  Constants 

=  3.8  x  102  lbfs/ft2 

R,„  =  53.3  ft  lbf/lbm  °R 

°R  =  °F  +  460 

1  cm  H20  =  2.045  psf 

1  psi  =  2.036  in  Hg 

q  =  AP/.93  (Calibration  factor) 
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Report  Requirements 

1.  Calculate  the  effective  Reynold's  [sic]  number  and  determine  the  corresponding 
empirical  drag  coefficient  from  Figure  22. 

2.  Utilizing  x-y  plotter  curves,  evaluate  the  integrals  of  equations  (4.1)  and  (4.2)  and 
determine  experimental  drag  coefficients  for  each  method,  applying  wake  and  solid 
blockage  correction  factors. 

3.  Determine  the  percentage  reduction  in  drag  afforded  by  a  turbulent  boundary  layer. 

4.  Briefly  discuss  results. 


Figure  22.  Cylinder  Drag  Coefficient 
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APPENDIX  B:  UPGRADED  LAB  HANDOUT 


CYLINDER  DRAG  EXPERIMENT 


Class  time:  :50 

Lab  time:  1:30 

Report  time:  3:00 

Total  time:  5:20 

♦♦♦Before  the  Lab*** 

You  will  need  to  read  this  handout  as  well  as  pages  11-15,  32-38  of  the  NPS  Laboratory 
Manual  for  Low  Sneed  Wind  Tunnel  Testing.  You  will  also  need  to  bring  either  3 'A " 
or  5‘4M  formatted  floppy  disk  with  at  least  100  KB  space  available. 

Objectives 

A  current  method  of  PC-based  automated  daia  acquisition  using  a  combination  of 
pressure  transducers,  potentiometers,  signal  conditioners  and  amplifiers  will  be  used  to 
collect  and  analyze  data  from  flow  over  a  non-rotating  cylinder.  Two  techniques  of  drag 
determination  (pressure  distribution  and  wake  analysis)  will  be  used,  demonstrating  basic 
wind-tunnel  methodology.  In  addition  to  these  two  primary  objectives,  students  will  be 
introduced  to  the  NPS  low  speed  wind  tunnel  operating  procedures  to  include  tunnel 
calibration  and  tunnel  correction  factor. 

Apparatus 

The  NPS  low  speed  wind  tunnel  is  described  in  detail  in  the  NPS  Laboratory  Manual  for 
Low  Speed  Wind  Tunnel  Testing.  Figure  21  shows  the  signal  path  for  the  data 
acquisition  system.  The  dimensions  of  the  cylinder  are  d=2.96  inches,  h=28.4  inches. 
The  dimensions  of  the  test  section  are  w=45.0  inches,  h(above  reflection  plate) =28.4 
inches.  There  is  a  single  pressure  port  near  the  center  of  the  cylinder.  The  cylinder  is 
mounted  to  an  electrically  rotatable  base  which  is  able  to  rotate  from  -20  degrees  to 
+  180  degrees  in  a  clockwise  fashion.  Behind  the  cylinder,  on  a  manual  traverse,  is  a 
combination  pitot-static  probe  apparatus.  The  probes  may  be  moved  along  the  traverse 
in  the  y-direction  from  0  to  +20  inches.  The  transducers,  potentiometers,  signal 
conditioners  and  amplifier  which  connect  these  measuring  devices  to  the  data  acquisition 
card  inside  the  PC  are  shown  in  Figure  21.  Specific  information  about  each  device  is 
located  in  Table  1.  The  wind  tunnel  will  be  run  in  4th  gear  at  15  cm  H20  Ap.  This 
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"Ap"  is  a  IS  cm  H20  difference  between  the  local  static  pressure  entering  the  contraction 
cone  (p,)  and  the  local  static  pressure  leaving  the  contraction  cone  and  entering  the  test 
section  (p2).  Previous  wind  tunnel  calibrations  have  experimentally  shown  this  particular 
difference  of  15  cm  H20  to  equate  to  approximately  115  MPH. 

In  general,  the  calibrated  NPS  low  speed  academic  wind  tunnel  equation  is 

Ap  -  0.243  (cm  Hfi)  *  0.895  0) 


Dynamic  pressure  is 


ifes*  section 


* 


(2) 


where 


AP  -  IP, -Pj(cm, [y_l(A)  f-JJL.  ,  ,3, 


and 


t  *  62.35  — 

'  water  ^ 


(4) 


Combining  equations  (1)  through  (4)  and  taking  density  in  slugs/ft3, 


! 

<1 

/-N 

T1 

2(Ap  -  0.243)  (S) 

0.895  p. 
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X>  TJ 


Terminology 

p0  =  total  (stagnation)  pressure 

p  =  static  pressure 

q  =  dynamic  pressure  =  p0  -  p 

(Po)*  =  freest  ream  total  pressure  (associated  with  freestream  velocity) 
=  freestream  static  pressure  (associated  with  freestream  velocity) 
=  freestream  dynamic  pressure  =  (po)„  -  p* 

(Po)w  =  wake  pressure  measured  behind  the  cylinder  at  the  probe 
((po)»  assumed  to  =  (po)*  when  taken  outside  wake) 
pw  =  static  pressure  measured  behind  the  cylinder  at  the  probe 
(pw  assumed  to  =  p„  when  taken  outside  wake) 
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qw  =  wake  dynamic  pressure  =  (p„)w  -  pw 
pcy)  =  pressure  measured  at  the  cylinder 
Dp  =  Form  Drag  (due  to  separated  flow) 

Df  =  Shear  Drag  (due  to  skin  friction) 

Cp  =  pressure  coefficient 

cdu  =  uncorrected  drag  coefficient 

cd  =  drag  coefficient  relating  to  total  drag  or  profile  drag 

cdp  =  drag  coefficient  due  to  Dp  only 

r  =  radius  of  cylinder 

9  =  azimuth  of  the  rotator  (signal  from  potentiometer) 

y  =  position  of  the  traverse  (signal  from  potentiometer) 

Conversion  Constants 

latitude  correction  =  -0.0245  in.  Hg 

temperature  correction  =  (see  Table  4-2  on  bulletin  board  near  tunnel) 
=  3.8  x  107  (lbf-sj/ft2 
R*r  =  53.3  (ft-lbf)/(lbm-°R) 

°R  =  °F  +  460 
1  ft  =  30.48  cm 
1  cm  H20  =  0.0142  psi 
1  psi  =  2.036  inches  Hg 

Qtest  section  =  (Ap  -  0.243)/. 895  (NPS  wind  tunnel  calibration  equation) 
TF  =  1.04  (NPS  wind  tunnel  turbulence  factor) 


TABLE  1.  Drag  Classification 


TOTAL  DRAG  -  D 

SHEAR  -  Df 
(due  to  skin  friction) 

Dominates  streamlined 
bodies 

NORMAL 
(due  to  pressure) 

Perm  - 

(due  to  $ej»aratk») 

Dominates  blunt  bodies 

Induced 

(due  to  lift) 

Parasite  Drag  =  Df  +  Dp  (3-D) 
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Pressure  Distribution  Method 

For  blunt  objects  in  subsonic  airspeeds,  form  drag  dominates  the  resultant  force  on  a 
non-rotating  cylinder  in  airflow.  The  pressure  distribution  method  only  measures  Dp. 
Using  Bertin  and  Smith’s  Aerodynamics  for  Engineers,  p.  60  or  Anderson’s 
Fundamentals  of  Aerodynamics,  p.  209: 


/•*  P~P  P~P 

c .  -  /  - -  cos0  60,  where  C„  =  - -  (6) 

dP  Jo  qm  P  qm 


Wake  Analysis  Method 

Comparatively,  the  wake  analysis  method  measures  not  force,  but  a  change  in  velocity. 
It  is  a  ratio  of  dynamic  pressures,  which  reduce  to  a  ratio  of  velocities  according  to 
equation  (2).  The  wake  analysis  method  accounts  for  profile  drag,  or  D{  +  Dp.  From 
Rae  &  Pope’s  Low-Speed  Wind  Tunnel  Testing,  p.  215: 


1  f  20 

C'-/o  ^ 

(Po)w  P w  (PqJw  p»  )  (7) 

Data  Acquisition 

The  heart  of  any  PC-based  data  acquisition  system  is  the  data  acquisition  A/D 
(Analog/Digital)  board.  This  device  takes  analog  signals  such  as  voltages  and  converts 
them  to  digital  signals  for  the  computer  software  to  use.  In  the  case  of  this  laboratory, 
pressure  signals  and  position  signals,  in  the  form  of  output  voltages  (Vdc)  are  first 
conditioned  for  optimum  signal  gain,  then  channeled  to  the  A/D  board  via  a  standard  RS- 
232  terminal.  The  A/D  board  simply  plugs  into  the  inside  of  the  PC  and  has  two  settings 
which  are  of  interest  to  the  user: 

1.  8/16  Channel  Select  -  8-channel  differential  input  (normal  or  floating)  or  16- 
channel  single-ended  input.  This  lab  uses  8-channel  floating  differential  input  (Channels 
0-7).  The  term  "floating"  merely  refers  to  the  method  of  providing  a  common  ground 
to  eliminate  unnecessary  noise. 


2.  Gain/ Range  -  manual  or  programmable.  The  A/D  board  used  in  this  lab  has  a 
programmable  gain/range  accessible  through  the  software.  The  gain/range  may  be  set 
differently  . for  each  input  device  and  for  all  uses  of  that  input  device.  Bipolar  and 
unipolar  refer  to  the  polarity  of  the  range.  An  example  of  bipolar  gain/range  would  be 
±5V.  Unipolar  might  be  0-5V. 

Once  the  signals  have  been  properly  fed  into  the  A/D  board  it’s  up  to  the  software 
to  use  the  signals  to  provide  actual  data  acquisition.  This  includes  accepting  the  data, 
displaying  the  data  and  storing  the  data.  Data  that  have  been  stored  (for  example  a  file 
of  two  columns  of  data:  elapsed  time  and  temperature)  may  then  be  retrieved  into  a 
spreadsheet  and  plotting  program  (such  as  Lotus  1-2-3*)  or  into  a  math  processor  (such 
as  MATLAB*)  to  be  analyzed.  The  convenience  of  the  PC  allows  the  operator  to  do  this 
at  one  place  and  at  essentially  one  time. 

A  typical  method  of  calling  a  data  file  into  MATLAB*,  for  instance,  would  be: 

>  > load  A:\pr*  dial. m 

>  >  whoa  %  will  display  list  of  matrix  |m.n]  (you  will  need  m) 

>  >  x  =  pre_disl(l:m.l);  %  create*  a  column  vector  x  of  the  tint  column  of  data  in  pre  dial 

>  >  y  -  prediald  %  create*  a  column  vector  y  of  the  second  column  of  data  in  pre  dial 

>  >  plot(x,y),grid  %  create*  a  plot  of  x  v*.  y 

In  this  AA2801  lab,  the  following  data  will  be  generated  by  LabTech  Notebook: 


TABLE  2.  File  Organization 


Column  1 

Column  2 

Column  3 

PreJDisl.m 

e 

Cp  (Raw) 

Cp  (theoretical) 

WakeAnl.m 

y 

> 

WakeAn2.m 

y 

> 

_  <Lw 
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If  the  A/D  board  is  the  "heart"  of  data  acquisition  then  the  software  is  the  "brains." 
Data  acquisition  software  programs  (such  as  LabTech  Notebook*)  allow  the  operator  to 
configure  "set-ups"  which  tell  the  PC  how  to  collect,  display  and  store  the  data.  Current 
software  uses  block  icons  which  may  be  cascaded  into  an  infinite  number  of  combinations 
of  configurations. 

Safety  Procedures 

Prior  to  executing  the  lab,  normal  pre-operating  safety  checks  should  be  made.  Students 
should  use  proper  hearing  protection  for  tunnel  speeds  above  5  cm  H20. 

Initial  Procedures  Prior  to  Running  the  Tunnel 

1.  Determine  atmospheric  pressure  (pItIB  =  (bare) + (temp.  corr.)+(lat.  corr.)) 

Pbaro  _  Patm  _ 


2.  Observe  tunnel  temperature...  T4lm  = _ 

3.  Ensure  the  following  prior  to  tunnel  runs: 

a.  Top  rack  signal  conditioners  -  ON  (switch/light  is  on) 

Caution 

Do  not  change  or  attempt  to  adjust  any  settings  at  this  time 

b.  Traverse  is  in  y  =  0.00"  position  (visual  inspection) 

c.  Rotator  is  in  8  =  0°  position,  toggle  switch  =  RUN,  speed  control  =  75 
(visual  inspection) 

d.  486DX/33MHz  PC  turned  ON  and  in  Windows  Program  Manager 

e.  HP  Laser  Jet  II  Plus  turned  ON  with  paper  loaded 

f.  Double  doors  and  door  to  room  H-033A  CLOSED 

4.  Double-click  "Build-Time"  under  the  "Wind  Tunnel"  menu  (you  are  now  in 
LabTech  Notebook’s  ICONview  -  you  may  need  to  click  the  uppermost  right  button 
to  stretch  the  window  so  that  you  can  see  the  entire  set-up) 

Note:  Because  of  the  memory  usage  involved,  waiting  for  WINDOWS 
to  complete  each  input  before  assigning  another 
will  reduce  the  potential  to  overload  the  PC;  if  the  PC 
is  overloaded,  chances  are  you  will  have  to  reset  and  start  over. 
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5.  Perform  ZERO  and  SPAN 


Calibration 

The  lab  technician  will  ensure  the  data  acquisition  lab  has  the  proper  initial  settings  and 
that  the  input  devices  are  properly  calibrated.  The  configuration  outputs  should  be 
"scaled"  at  this  time  using  the  set-up  "ZERO"  and  "SPAN." 
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TABLE  3.  AA2801  Cylinder  Lab  Initial*  Settings 


Channel 

0 

1 

2 

3 

4 

Run 

1 

2,3 

2,3 

1 

2,3 

Signal 

PcyrP  oo 

(Po)w'Pw 

y 

e 

(Po)oo‘Poo 

HI 

37 

36 

35 

34 

33 

LO 

18 

17 

16 

15 

14 

LLGND 

19 

28 

29 

15  -J->  29 

19 

Sig.  Cond, 

#2 

#1 

#3 

Pac.  #2 

m 

Amp? 

No 

No 

No 

Pac.  n, 
xlOO 

No 

Zero  pt. 

0  cm 

0  cm 

y=0" 

II 

I 

•— * 

OO 

o 

0  cm 

‘Zero  (v) 

O.OOv 

O.OOv 

O.OOv 

-0.1 8v 

O.OOv 

Span  pt. 

15  cm, 
0=0°,  y=0" 

15  cm, 
y=0" 

o 

rs 

II 

>, 

0=  +  18O° 

15  cm 

‘Span  (v) 

0.009v 

0.0054v 

.0200v 

1.80v 

0.009v 

Scale  Factor 

xl  10 

xl85 

xlOOO 
+  2.0 

xlOO 

xl  10 

Polarity 

Bipolar 

Bipolar 

Bipolar 

Bipolar 

Bipolar 

Gain 

±0.625v 

±0.625v 

±0.625v 

±2.50v 

±0.625v 

"12  Bit  Resolution” 

(2n  =  4096  divisions  available) 

example: 

vo  —  =  ±  0.305  mV  resolution 

(4096  -  1) 
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*Table  3  show  "initial  settings"  only  and  is  intended  for  the  technician  to  use  prior  to  the 
beginning  of  the  AA2801  Cylinder  Lab.  LabTech  Notebook  (LTNB)  setups  "ZERO," 
"SPAN"  and  accompanying  directions  should  be  used  for  final  settings. 

Using  '’ZERO”  and  "SPAN” 


All  input  devices  are  referenced  to  some  zero  and  to  some  secondary  value,  which  yields 
a  "span."  That  is,  some  mode  of  their  output  where  the  user  wants  0.00  volts  and  some 
other  voltage  as  a  return.  Everything  is  considered  linear  after  the  zero  and  span  are  set 
to  a  desired  output.  Normally  a  Digital  Volt  Meter  (DVM)  is  used  to  verify  such  an 
output.  However,  for  automatic  data  acquisition,  it’s  more  important  to  verify  the  output 
at  the  PC  since  the  date  *re  being  collected  here.  (It  is  possible  for  noise  to  affect  what 
is  read  after  the  signal  conditioner  so  that  the  signal  has  changed  by  the  time  it  is 
displayed  on  the  monitor.)  The  "ZERO"  and  "SPAN"  set-ups  were  designed  to  give  the 
user  a  visual  display  of  his/her  zero  and  span  for  each  device  used  in  the  cylinder  lab. 


Figure  22.  ZERO  Showing  Signal  Figure  23.  ZERO  Showing  Signal  After 

Somewhere  Less  Than  Zero  Volts  Being  "Flown"  Up  to  Zero  Volts 

Here  are  the  steps: 

1 .  Using  the  486  PC  on  the  right,  RECALL  ZERO  in  LabTech  Notebook  Build-Time. 

2.  Ensure  the  tunnel  is  not  running  by  verifying  the  tunnel  water  micromanometer 
does  in  fact  read  0.00  cm  with  the  bottom  of  the  water  meniscus  just  resting  on  the 
cross-hair.  Verify  at  the  micromanometer,  not  the  PC. 

3.  Select  RUN. 

4.  The  technician  has  already  pre-configured  the  signal  conditioners  to  yield  a  ballpark 
zero.  Using  the  display  presented,  unlock  and  turn  each  "zero"  knob  on  the  upper 
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signal  conditioners  (#1,2  and  10)  until  the  colored  "meatball"  is  centered.  This 
is  very  easy  for  Naval  aviators!  A  line  display  reflects  the  actual  zero  line  for  each 
device.  Because  of  the  complexity,  the  technician  will  ensure  the  rotator  and 
traverse  are  properly  scaled. 

5.  To  illustrate  the  effect  of  undesired  noise  on  a  data  acquisition  system,  turn  the 
DVM  on  and  plug  the  red  probe  into  the  green  outlet  of  the  second  signal 
conditioner;  plug  the  black  probe  into  the  white  outlet.  Wait  for  30  seconds.  What 
happens  to  your  display?  (Don’t  forget  to  unplug  the  DVM  probes!) 


6.  Exit  ZERO  and  repeat  step  1  for  SPAN. 

7.  After  the  operation  brief  and  safety  check  of  wind  tunnel,  run  the  tunnel  up  to  15 
cm  H20. 

8.  Select  RUN.  This  time,  the  signal  will  be  referenced  to  1.0  (unity). 

9.  Again,  the  technician  has  already  pre-configured  the  signal  conditioners  to  give 
ballpark  voltage  output  which,  when  scaled  by  the  set-ups  in  LabTech  Notebook, 
will  yield  1.0.  Using  the  display  presented,  unlock  and  turn  each  "span”  knob  on 
the  upper  signal  conditioners  (#1,2  and  10)  until  the  colored  "meatball"  is  again 
centered.  Another  line  display  reflects  the  actual  1 .0  line  for  each  device.  Again, 
the  technician  only  will  ensure  the  rotator  and  traverse  are  properly  scaled. 

10.  Press  "REV”  on  the  rotator  control  box  to  reset  the  cylinder  to  0  =  -20°. 

1 1 .  The  setups  are  now  properly  scaled.  Leave  the  tunnel  running  and  proceed  to  step 
1  of  the  cylinder  lab. 

Lab  Procedures 

RUN  #1 

1.  RECALL  PRE  DISl 

2.  Single-click  the  "Run"  menu  and  single-click  "Run” 

(you  will  see  a  window  which  asks  you  if  you  want  to  "Overwrite  previous 
log(s)?") 

3.  Single-click  "Yes"  and  expand  the  window  (the  display  configured  for  this  third  of 
the  lab  is  now  present  -  data  acquisition  has  been  "triggered"  to  begin  when  the 
input  from  the  azimuth  block  passes  through  -17°  for  the  first  time) 

4.  Momentarily  push  the  black  "FWD"  button  on  the  rotator  control  box,  observe  the 
following: 

a.  the  cylinder  rotates  from  -20  to  +180  degrees 

b.  as  the  cylinder  passes  -17  degrees,  data  is  displayed  on  the  PC: 

Azimuth  and  current  time  displayed  in  meter  boxes 
X-Y  plot  showing 
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Green  line:  theoretical  Cp  vs.  6 

Cyan  line:  actual  raw  Cp  vs.  6  (collected  at  Sampling  Rate  of  10  Hz) 
Black  line:  actual  smoothed  Cp  vs.  9  (moving  average  with  parameter 
of  10  points) 

c.  boundary  layer  separation  occurs  (the  wake  begins) 

d.  data  acquisition  stops  as  cylinder  passes  through  +180  degrees  but  the 
display  remains  on 

5.  To  make  a  print-out  of  the  display  hit  ALT-Print  Screen  (this  sends  a  BitMap  to 
the  ClipBoard  in  Windows*) 

6.  Single-click  "File”  then  single-click  "Exit"  to  leave  the  display  and  return  to 
ICONview  (you  will  be  asked  "Do  you  really  want  to  exit?"  -  single-click  "OK") 

7.  Exit  ICONview  by  single-clicking  "File,"  "Exit,"  and  "OK"  in  succession 

8.  Double-click  "paintbrush"  in  the  Windows*  Accessories  menu 

9.  Single-click  "View,"  then  deselect  (by  single-clicking)  "Tools  and  Linesize,"  and 
"Palette" 

10.  Click  the  uppermost  right  button  to  stretch  the  paintbrush  window  to  its  limits 

1 1 .  Select  "Edit"  and  "Paste"  (a  copy  of  the  previous  display  will  appear) 

12.  Select  "File"  and  "Print" 

13.  Reduce  the  scaling  to  80%  and  hit  "enter"  (the  printer  will  now  print  a  copy  of  the 
display) 

14.  Select  the  uppermost  "File"  and  "Exit"  (paintbrush  will  ask  you  if  "you  want  to 
save  current  changes?  -  single-click  "No") 

RUN  n 

15.  Repeat  steps  1,2  substituting  WAKE  AN1  for  the  first  set-up 

16.  Single-click  "Yes"  and  expand  the  window  (the  display  configured  for  this  third  of 
the  lab  is  now  present  -  data  acquisition  has  been  "triggered"  to  begin  when  ALT-1 
is  pressed  on  the  keyboard. .  .press  ALT-1  when  the  next  step  is  started) 

17.  Have  someone  begin  cranking  the  traverse  handle  at  a  quick,  even,  steady  rate  and 
observe  the  following: 

a.  the  traverse  moves  from  0  to  20  inches 

b.  after  ALT-1  is  pressed,  data  begin  to  be  displayed  on  the  PC: 

Position  and  time  are  displayed  in  meter  boxes 

X-Y  plot  showing 
Blue  line:  smoothed 


vs.  y 
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(moving  average  with  parameter  of  25  points  and 
collected  at  Sampling  Rate  of  10  Hz) 

c.  the  position  and  relative  strength  of  the  wake 

d.  data  acquisition  continues  as  traverse  passes  through  19  inches  and  display 
remains  on 

18.  Press  ALT-Print  Screen  to  save  display,  then  select  File  and  Exit  immediately  to 
stop  data  acquisition.  To  make  a  print-out  of  the  display  repeat  steps  5-14 

19.  Stop  the  wind  tunnel  at  this  time  using  tunnel  operating  procedures,  open  test 
section 

RUN  #3 

20.  Repeat  ZERO  to  recalibrate  your  signals,  if  necessary. 

21 .  Apply  two  1/2"  x  20"  x  1  mil  (approx.)  pieces  of  cellophane  tape  to  cylinder  1  % 
inches  on  either  side  of  the  static  port  (this  will  trip  the  boundary  layer,  creating 
turbulence  and  thus  delaying  boundary  layer  separation) 

22.  Crank  the  traverse  back  to  y  =  0  and  close  up  the  test  section 

23.  Run  the  tunnel  up  to  15  cm  H20  again 

24.  Repeat  SPAN  to  recalibrate  your  signals,  if  necessary. 

25.  Repeat  steps  15-19  substituting  WAKEAN2  for  the  second  setup  (notice  that  the 
wake  size  is  reduced) 

26.  Note  the  tunnel  temperature  at  shut  down:  T  = _ 

27.  You  are  complete  -  Exit  LabTech  Notebook  and  double-click  "Exit  Windows"; 

28.  Place  a  formatted  disk  in  either  floppy  drive.  If  you  are  using  a  3 Vi"  disk,  at  the 
DOS  prompt  type  (otherwise  substitute  A  for  B): 

C:\  >  copy  C:\windtunl\pre_disl.m  B:\predisl.m  (ENTER) 

C:\>  copy  C:\windtunl\wake_anl.m  B:\wake_anl.m  (ENTER) 

C:\>  copy  C:\windtunl\wake_an2.m  B:\wake_an2.m  (ENTER) 

This  will  copy  the  three  data  files  onto  your  disk  for  later  analysis. 

29.  Secure  the  CPU,  Monitor  and  Printer;  ensure  the  wind  tunnel  is  properly  secured. 


Report  Requirements 

1 .  Using  a  data  analysis  program  such  as  MATLAB*  and  data  from  the  lab,  calculate 
the  integrals  of  equations  (1)  and  (2).  Be  sure  to  apply  proper  wake  and  solid 
blockage  correction  factors  to  find  experimental  drag  coefficients  for  each  of  the 
three  runs.  See  pages  35-38  of  the  NPS  Laboratory  Manual  for  Low  Speed  Wind 
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Tunnel  Testing  for  more  information.  Explain  in  detail  how  you  found  the  drag 
coefficients.  List  any  equations  used  and  how  you  used  them. 

2.  Using  the  cdu  from  Run  #1,  calculate  the  effective  Reynolds  number  (Re^p)  and 
determine  the  corresponding  empirical  drag  coefficient  using  Figure  24  (from 
Anderson’s  Fundamentals  of  Aerodynamics,  p.  229)  and 

-  KtmmmJLl+t)-7F  <*> 


3.  Briefly  discuss  the  differences  between  the  pressure  distribution  results,  the  wake 
analysis  results  and  the  empirical  drag  coefficient. 

4.  Determine  the  percentage  reduction  in  drag  afforded  by  an  induced  turbulent 
boundary  layer.  List  at  least  two  real-life  examples  of  how  turbulence  is  used  to 
delay  the  onset  of  drag  (boundary-layer  separation). 


81 


82 


85 


87 


88 


89 


90 


APPENDIX  D:  MATLAB  CODE  FOR  RETRIEVAL/DATA  ANALYSIS 
Filename:  DRAGl.m 


%  DRAG1.M  needs  3  files  on  s  3-1/2"  disk  labeled: 

%  B:\pre  disl.m;  B:\wske  sal. m;  B:\wakc_an2.m 
%  If  you  haven't  placed  a  disk  in  the  B-drive  with" (best  three  files,  as 
%  named,  do  so  at  this  time.  If  you  have  your  files  on  a  3-1/4"  disk, 

%  you  will  need  to  get  into  this  program  "DRAGl.m"  and  edit  the  lines 
%  with  the  load  cosmnsods  to  read  A  vice  B,  for  A-drive. 

%  Remember,  the  student  version  of  MATLAB  will  not  execute  s  file  with 
%  more  than  1024  elements  The  Wake  An  files  have  at  least  that  many. 

%  Next,  it  loads  these  three  files  into  MATLAB  and  then  begins  to 
%  break  them  down  into  column  vectors  for  use  in  a  modified  Trapezoidal 
%  rule.  Instead  of  having  equiapaccd  x  points,  as  required  by  the 
%  trapezoidal  rule,  the  difference  between  each  successive  set  of  x-values 
%  are  multiplied  by  the  average  between  each  associated  set  of  y-values. 

%  These  products  are  then  summed.  All  cd’s  are  then  corrected  for  blockage. 

%  The  following  abbreviations  are  used: 

%  cdpdu  =  s  uncorrected  cd  using  eqn  2  (for  use  in  Re(cff)  calculation) 

%  Reeff  -  =  effective  Reynolds  number  for  use  w/  Anderson's  Figure  3.39 
%  cdpdl  =  =  cd  from  pressure  distribution  method;  eqn  I  -  equiapaccd  assumed 
%  cdpd2  =  =  cd  from  pressure  distribution  method;  eqn  2  -  modified  trap,  rothd 
%  pftl  =  =  a  poly  fit  of  RUN  01  values  was  done  and  a  eqn  1  was  applied 
%  cdwal  =  =  drag  coefficient  (RUN  #2)  via  the  wake  analysis  method;  eqn  2 
%  cdwa2  =  =  drag  coefficient  (RUN  03)  via  the  wake  analytia  method;  cqn  2 
%  sepaz  =  =  the  azimuth  in  degrees  where  separation  occurs  (estimate) 

%  wakebeg  =  =  the  estimated  beginning  of  the  wake  (in  y) 

%  wakccnd  =  =  the  estimated  coding  of  the  wake  (m  y) 

%  percrcdu  -  —  the  %  drag  reduction  afforded  by  tripping  the  boundary  layer 
%  Expect  some  plots  after  the  program  finds  each  cd.  Hit  any  key  to  continue 
%  after  the  computer  pauaes  with  each  plot  or  cd.  At  this  time  type: 

%  dragl  after  the  >  >  prompt.  Watch  the  B-dnve  light  to 

%  clue  you  as  to  the  proximity  to  your  first  value.  It  takes  shout  73  seconds. 

load  B:\prc  disl.m 

load  B:\wake_anl.m 

load  B:\wske  an2.m 

|ml,nl|-size(pre  dial); 

j  m2  ,n2 1 = sizc(wa£e_an  1 ) ; 

j  m3  ,n3  j  =  size(wakc~an2) ; 

xl  =  prc_disl(l:ml,T); 

yl  =pre  disl(l:ml,2); 

x2  =  wa£e  anl(l:m2,l); 

y2~wakc~anl(l:m2,2); 

x3  =wake_an2(l  :m3,l); 

y3  =  wake  an2(l:m3,2); 

xxl  =xl.*fpi/180);  *  degree*  to  radian  conversion 

yyl  =yl.*cos(xxl);  %  each  elemental  Cp  is  multiplied  by  its  associated 

%  cos(thcta) 

%  approx,  cquispace  by  letting  delta  x  =  total  length/#  samples;  then 

*  use  trapezoidal  rule  (Holman,  eqn.  3-16,  p.  193);  I  call  this  "eqn.  1  ": 
cdpdl  =((yyl(ml)-yyl(l))/2-i-(sum(yyl)-(yyl(ml))))*(xxl(ml)/(ml-l)); 

%  this  is  only  slightly  different  from  eqn.  1;  here  I  do  not  rely  on 

%  equiapaccd  x-values;  instead  1  sum  each  successive  product  of  dx  and  dy; 

%  I  call  this  "cqn.  2": 
for  k=2:ml; 

s0(k) =(yyl(k')+yyl(k))/2  *abs(xx  1  (k)-xx  1  <k- 1 )) ; 

end 

cdpd2 =sum(sQ); 

pi  = poly  fit(  xxl  ,yy  1 ,5);  %  this  makes  a  polyfit  of  the  values  given  to  a  3th 

*  degree  (order)  polynomial 

f=po!yvai(pl,xxl);  %  this  evaluates  that  polynomial  at  each  x-value 
pftl  =((f(ml)-f(l))/2+(sum(f)-(f(ml))))*(xxl(ml)/(ml-l));  %  eqn  1 
for  i=2:m2; 

sl(i)=(y2(i-l)+y2(i))/2*abs(x2(i)-x2(i-l));  %  eqn  2 

end 

cdwal  =sum(sl)*2/2.96; 
for  j=2:m3; 

s2(j) = (y3(j- 1 )  + y  3(j))/2  *abs(x3(j)-x3(j- 1 »;  *  eqn  2 

end 

cdwa2  =  sum(s2)  *2/2 .96; 
cdpdu =cdpd2; 
cdpdu 

Reeff =270477 .9*1. 04*(1  +(. 00234 +  .0329*cdpdu)); 

Reeff 

%  this  next  calculation  adds  the  wake  and  solid  blockage  factors 
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cdpdi=cdpdlNI-3*«0.J2M9J.43)/(S.*M44)*l.J>-2*<(2.96«cdpdl)/<2«43)»; 

cdpdl 

cdpd2=cdpd2*(l-3»«0.J2*I93.43>/(».l»l44)*l.J)-2*t<2.96«cdpd2)/(2MJ))>; 

cdpd2 

pftl-pftl*(l-3*«0.J2M9J.43)/<*.»M44)Al,J)-2*((2.96»pftl)/<2*4J)»; 

pftl 

pAUK 

plot{xxl,yyl,xxl,f),fnd,Mle('Cp*coe<thcta)  vs  theta'),xlebeK'theU')4tex«'Jth  order  polynomial') 
pause 

cdwal  “  cdwal  •(l-3*((0.32,l93.43)/(E.E*l44)AI.J)-2*((2.96*cdwal)/(2*4S))); 
cdwal 

%  what  ia  the  anumali  (location)  where  separatioa  occurs? 
for  t>l:ml; 
if  yl(t+3)  >  yl(t); 

sep  =»  xx  1  (t) ; break  .cod 
end 

yp=2*2.9h*sin(sep); 
tepaz=tcp*l  SO/pi; 
sepaz 

plo«xl,yl,scpaz,yl(l),’,’),*nd.title(’Cp  vs  theta ’),xiabcl('tbeta'),xtext('ieparatica') 
pause 

wakebec  “  10-(yp/2); 

wakebcc 

wakecnd-10+(yp/2); 

wakeend 

pause 

plot(x2,y2),(nd,titie('ij  ratios  diff  vs  y'),xUbel('y') 
pause 

cdwa2=cdwa2»<l-3*«0.52»195.43)/<8.«»144)Al.J)-2*«2.96»cdwa2)/(2»45))); 

cdwa2 

percrcdu  =  I00-(cdwa2/cdwal)*100; 

pcrcredu 

pause 

ploKx3,y3,x2,y2),(rid,tidc('a  companion'), xlabel('y’),gicxl(Tnppcd  Boundary  Layer’) 

pause 

end 
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APPENDIX  E:  MATLAB  CODES  FOR  TUNNEL  VELOCITY  TABLES 
Filename:  SPuEDS.m 

%  Thu  Mallab  file  SPEEDS. m  crease*  conversion  column*  in  different 
%  uni's  lor  the  NPS  Low  Speed  Wind  Tunnel;  (he  equations  herein 
%  ere  baaed  on  data  obtained  durum  calibrations  made  by  Professor 
%  R  Howard  and  LCDR  Clay  MiUer  on  26  April  1993. 

■4 

%  Thu  program  is  written  out  to  show  conversions  m  approximately 
%  I  page  formats  utd  will  pause  and  wait  for  any  keystroke 
%  between  displays  to  allow  user  to  print  data 
% 

%  You  will  need  to  remember  the  following  order  of  columns 
* 

%  DEL  P  ft/ see  m/s  ku  mph  km/hr 
* 

1=  IS:. 03:1.0;  %  Thu  is  the  ref  cm  of  H20  on  the  nucromanomeur 
%  in  ft/scc: 

f=sqri((2*(i.  *(62. 35/30.3)-. 24326))/(0. 89467*. 002378)); 

%  in  m/s: 

m=sqrt(  (2  *(i.  *(62. 35/30.5)-. 24326))/(0. 89467*. 002378))*.  3048; 

*  “*  k  =  «qrt((2*(i.  *<62. 35/30.5)-. 24326))/(0.S9467*. 002378))*. 5921 ; 

%  in  mph: 

mi =sqit((2*(i.  *(62. 35/30.5)-. 24326))/<0. 89467*. 002378))*. 681 8; 

%  in  km/hr: 

km“sqiK(2*(i.*(62.35/30.5)-.24326))/<0.89467*.002378))*l  .097; 
s  =  (i'  f  m'  k’  mi'  km'l 


plot(i,f).grid.xlabeI('H20  cm'),ylaheI('ft/sec').titlr('velociiy') 
pause 

plot(i,m),gnd,xlabcl(’H20  cm'),ylabel('m/s'),titk('vclociiy') 


pka(i,k),gnd.slsbel(,H20  cm’),yiabcK’ku’).litlc4'vclocity') 


plot(i.mi).gnd,xJabcl(’H20  cm').ylabelCinpli').title('velocity') 
pause 

plot(i.km),grid.xlabcl('H20  cm'),ylabel('kin/hr'),title('velocity') 
pause 
end 
% 

i  =  I  0:  5:  10.0;  %  This  is  the  ref.:  cm  of  H20  on  the  mkromanomclcr 
%  in  ft/scc: 

f =sqrt(  (2  *(i.  *(62-35/30.3)-. 24326))/(0. 89467*. 002378)); 

%  in  m/s: 

m =sqrt((2  *(i.  *(62.33/30.3)-. 24326))/(0.89467*. 002378))*. 3048; 
%  in  kU: 

k*aqrt((2*(i.  *(62. 33/30.3). 24326))/(0.89467*. 002378))*. 5921; 
%  in  mph: 

mi =sqrt((2*(i.  *(62. 35/30.5)-. 24326))/(0. 89467*. 002378))*. 681 8; 

%  in  km/hr: 

km=sqrt((2*(i.  *(62. 35/30.5)-. 24326))/(0.89467*.002378))*l  .097; 
s  =  (i'  P  m'  k'  mi'  km'l 
pause 

pkH(i,D,grid,xlabel('H20  cm’),ylabel(’ft/ sec'), Ullcl ’velocity’) 
pause 

plot(i,m).gnd,xlabel('H20  cm’),yixbel('m/t’),Utie('velocity’) 
piuac 

plot(i,k).grid.xlabel('H20  cm').ylabel(‘kU’),tiUe(‘velocity’) 
pause 

plot(i,mi),gnd,xlabcl('H20  cm'),ylabel('mph’),titk('veloeity') 
pause 

plot(i,km),grid,xlabel(’H20  cm'),ylabel('km/hr’),tiae('velocity') 
pause 
end 
% 

i*10.0:. 3:20.0;  %  This  is  the  ref. :  cm  of  H20  on  the  micromanometer 
%  in  ft/sec: 

f=  n|it((2*(i.  *(62. 35/30.5)-. 24326))/(0.89467*. 002378)); 

%  in  m/s: 

m  =  sqtt((2*(i.  *(62.35/30.5)-. 24326))/(0.89467*.002378))*.  3048; 

m  k=sqrtf{2*(i. *(62.35/30.5)-.24326))/(0.89467*.002378))*  5921 ; 
%  in  mph: 
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mi  =  sqrt((2*(i.*(62.35/30.5)  .24326))/(0.*9467*.00237*))*  6*1*, 
%  in  km/fer: 

km  =  miK(2*(i.  *(62. 33/30.3)-. 24326))/(0.  *9467*. 00237*))*  1.097; 
s=|i'  P  m  k'  mi*  km‘| 

pmt 

pk*(i,0,fnd,ilabei(‘H2O  cm').ylabel('ft/sec’),utk('velocity’) 

pause 

plosion)  ,srid,xiabcl('H30  cm’), ylnbcK  W»'),tnk(  velocity') 
pause 

p)ot<i,k),Xnd,iUbclCH20  cm'),yiabel{'kti'),UUc(ve)ocity’) 

pause 

pka(i,mj),|nd,»Ubcl('H20  cm'), ylabcJ('mpb'),titlc< 'velocity') 

pime 

pk*(i,km).{nd,]Uabcl('H20  cm'),ylabel<'kjn/lu’),utk('velocity') 
pause 
end 
* 

i  =  20.0:. 5:30.0;  *  Thii  u  the  ref.:  cm  of  H20  on  the  fiucromanoraeUr 
%  in  fi/sec: 

f=*lit((2*(i.  *(62. 35/30.5)-.  24326))/(0.  *9467*. 00237*)); 

%  in  m/i 

m  =  M|it<(2*(i.  *(62. 35/30.5)-.  24326))/(0. 19467*. 002371))*.  304t; 
%  in  kti: 

k  =  «qrt((2*(i.  *(62.35/30.5)-. 24326))/(0.*9467*.00237*))*. 5921; 

%  m  mph: 

uu=sqil((2*(i.*(62  35/30.5)-. 24326))/(0. *9467*. 00237*))*. 6*1  S; 

%  in  km/hr: 

km  =  sqit(  (2  *(i.  *(62. 35/30.5)-.  24326))/(0.  *9467*. 00237l))*l  .097; 
i  =  |i'  f  m'  k'  mi*  km') 
pause 

pk4(i,f),(nd.xlabel/'H20  cm'),ylabel('ft/tec'),titie< 'velocity') 
pause 

plot(i,m),fnd,iUabel(’H20  cm'),ylabel(‘m/t'),title(‘velocity') 
piuae 

ploM i ,k),gnd,xUbcl< * H20  cm  *)  ,y  labcl< ’kti  *)  ,UUc< ’veloc aty * ) 
pause 

pkH( i ,mi ) ,gnd,  xJabeK * H20  cm’),yUbel(,Inph’),Utk(,veloc^ty,) 
pause 

pk*(i,kjn),gnd,xlabcl<*H20  cm'),yUbcl(’km/hr’),uUc(’vclocity’) 

pause 

end 
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Filename:  USER  SPD.m 

%  Thi»  Mallab  file  USER  SPD.m  allows  a  user  to  provide  MAT  LAB  with 
%  a  desired  velocity  and  Save  MATLAB  return  a  setting  for  the  water 
%  nucroma&oavctcr  (in  cm  H20)  to  be  used  with  the  KPS  Low  Speed  Wind 
%  Tunnel;  the  equations  herein  are  based  on  data  obtained  during 
%  calibrations  made  by  Professor  R.  Howard  and  LCDR  Clay  Miller  on 
%  26  April  1993. 

* 

% 

i=mput('MPH?  y/n  |n|:\’s’); 
if  i  =  =  'y' 

k  “inputl’ Enter  speed  desired  with  a  decimal  point  ’) 
DEL_P=((k./.6818)*2*(.89467*.002378/2)  +  .24326)*(30.S/62.35); 
('You  need  to  set'), DEL  P 
pause 

j=0.5:.5:40.0; 

p  =  sqrt((2*(j  *62. 35/30.5-0.24326))/(. 895*. 002378))*. 6818; 
plot(j,p,DEL  P,k,'*'),grid,xlabel('H20  cm’),ylabel(’MPH‘),title(’velocity’) 
pause 
break  .end 
elscif  i  =  =  'n' 

i = input!' Ft/sec?  y/n  |n|:',’s’); 
if  i=  =  ’y' 

k=input('Enter  speed  desired  with  a  decimal  point:  ') 
DEL_P=(k.A2*(.89467*.002378/2)  +  .24326)*(30.5/62.35); 

(’You  need  to  set'), DEL  P 
pause 

1=0.5:5:40.0; 

p=sqrt((2*(j*62.35/3O.5-0.24326))/(. 895*. 002378)); 
plot(j,p,DEL  P,k,'*’),gnd,xlabcl(’H20  cm’),ylabel('Ft/sec'),title(’velocity') 
pause 
break.end 
elseif  i=  =  V 

i  =  input(’m7s?  y/n  (nl:',’s'); 
if  i=  =  'y‘ 

k  =  input( 'Enter  speed  desired  with  a  decimal  point:  ') 

DEL  P=((k./.3048)*2 *(.89467*. 002378/2) +  .24326)*(30.5/62. 35); 

('You  need  to  set 7.DEL  P 

pause 

j  =  0.5:. 5:40.0; 

p = sqrt((2*(j*62  35/30.5-O.24326))/(.  895*.002378))*.3048; 
plot(j,p,DEL_PIk,’*').grid,xlabel('H20  em’),ylabel(’m/s’),title(’velocity’) 
pause 
break.end 
elseif  i=  =  'n’ 
i  =  input(’kts?  y/n  lnl:',’s'); 
if  i  =  =  'y ' 

k  =  inputC Enter  speed  desired  with  a  decanal  point:  ’) 

DEL_P=((k./.5921)A2*(. 89467*. 002378/2) +  .24326)*(30.5/62. 35); 

(’You  need  to  set'),DEL_P 

pause 

j=0.5:. 5:40.0; 

p=sqit((2*(i*62.35/30.5-0.24326»/(. 895*. 002378))*. 5921; 
plot(j,p,DEL  P,k,'*'),grid,xlabel(’H20  cm'),ylabcl('kts’),title('velocity’) 
pause 
break.end 
elseif  i  =  =  V 

i  =  inputCkm/hr?  y/n  |n|:’,’s'); 
if  i  =  =  'n' 

k  =  inputC  Enter  speed  desired  with  a  decimal  point:  ') 

DEL_P=((k. /I. (w7)',2*(. 89467*. 002378/2) +  .24326)*(30. 5/62. 35); 

('You  need  to  set'),DEL  P 

pause 

{=0.5:. 5:40.0; 

p=sqrt((2*(j*62.35/30.5-0.24326))/(.895*.002378))*l  .097; 
plot(j,p,DEL  P.k,'*'),gn<t.xlabel(’H20  cm'),ylabel(’km/hr'),title('velocity') 
pause 
break.end 
else  i=  =  ’n' 
end 
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